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ABSTRACT 


The  Marine  Corps  uses  its  Selective  Reenlisted  Bonus  (SRB)  Program  to 
influence  Marines  to  reenlist  for  a  designated  tenn  into  certain  Military  Occupational 
Specialties  (MOSs)  in  order  to  reach  planned  manpower  goals.  The  bonus  amount  is 
determined  by  selecting  an  “SRB  multiplier”  for  each  combination  of  MOS  and  Zone 
(“MOSZ”).  (“Zone”  corresponds  to  length  of  service.)  A  higher  multiplier  means  a 
larger  bonus  and  leads  to  a  higher  percentage  of  Marines  reenlisting.  That  percentage, 
predicted  by  an  existing  forecasting  model,  is  assumed  exact  here. 

The  “SRB  multiplier  model”  assigns  multipliers  to  minimize  a  sum  of  weighted 
squared  deviations  from  MOSZ  targets  subject  to  a  budget  constraint.  This  model  is 
implemented  as  a  generalized  assignment  problem,  and  solved  approximately  on  a 
personal  computer  using  Lagrangian  relaxation  and  a  secondary  heuristic.  (The 
algorithm  is  programmed  in  Visual  Basic  for  Applications  and  has  an  Excel  interface.) 

Data  for  FY04  shows  491  bonus-eligible  MOSZs.  With  up  to  11  possible 
multiplier  values,  this  yields  a  model  with  5,401  0-1  variables  and  491  constraints.  A 
solution  within  0.0018%  of  optimality  is  reached  in  1.4  seconds  on  1.58  GHz  personal 
computer.  Standard  integer-programming  software  verifies  the  correctness  of  the 
solution. 
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EXECUTIVE  SUMMARY 


The  United  States  Marine  Corps  (USMC)  uses  its  Selective  Reenlisted  Bonus 
(SRB)  Program  to  influence  Marines  to  reenlist  for  a  designated  term  into  certain 
Military  Occupational  Specialties  (MOS)  in  order  to  reach  planned  manpower  goals.  The 
amount  of  the  bonus  is  detennined  by  selecting  an  SRB  multiplier  for  each  combination 
of  MOS  and  Zone  (“MOSZ”).  (“Zone”  partitions  Marines  into  three  groups  based  on 
current  length  of  service.)  A  higher  multiplier  means  a  larger  bonus  and  leads  to  a  higher 
percentage  of  eligible  Marines  reenlisting.  That  percentage  is  predicted  by  an  existing 
forecasting  model,  whose  predictions  are  assumed  exact. 

The  total  amount  of  SRB  funding  available  for  a  fiscal  year  (FY)  is  limited  by  the 
amount  approved  and  allocated  in  the  FY  Defense  Budget.  Thus,  not  all  MOSZs  can  be 
assigned  the  multiplier  necessary  to  reach  reenlistment  goals.  Hence,  it  is  important  to 
allocate  the  available  SRB  funding  to  those  MOSZs  that  are  critical  to  accomplishing  the 
Marine  Corps’  mission.  So,  the  problem  that  is  addressed  by  this  thesis  is  to  detennine 
the  optimal  set  of  multipliers  to  reach  reenlistment  goals. 

The  “SRB  multiplier  model”  assigns  multipliers  to  minimize  a  sum  of  weighted 
squared  deviations  from  MOSZ  targets  subject  to  a  budget  constraint.  This  was  described 
and  implemented  as  a  generalized  assignment  problem  in  a  1986  Naval  Postgraduate 
School  thesis  by  U.S.  Marine  Corps  Captain  Dean  D.  DeWolfe,  entitled  Determination  of 
Selective  Reenlistment  Bonus  Multipliers  in  the  United  States  Marine  Corps.  This  model 
has  not  been  used  by  the  USMC,  however,  because  it  is  no  longer  current  with  respect  to 
SRB  rules  and  parameter  definitions,  and  because  it  was  not  implemented  in  a  user- 
friendly  programming  environment. 

This  thesis  overcomes  the  earlier  problems  with  the  SRB  multiplier  model  by  (a) 
updating  parameter  definitions  and  other  model  constructs  to  represent  the  current  SRB 
environment,  and  (b)  creating  a  flexible  implementation  of  DeWolfe’s  Lagrangian-based 
solution  algorithm  on  a  personal  computer.  The  flexible  implementation  should  allow  the 
model  and  solution  algorithm  to  track  changes  in  data  and  structure  for  the  SRB  program 
as  it  evolves  in  the  future.  The  algorithm  is  programmed  in  Visual  Basic  for  Applications 


xv 


and  is  supplied  with  a  user  interface  through  Microsoft  Excel.  The  algorithm  first 
optimizes  a  Lagrangian  function  that  relaxes  the  budget  constraint,  and  saves  the  best 
feasible  solution  it  finds  during  this  process.  The  algorithm  then  attempts  to  improve  that 
solution  with  a  greedy  marginal-rate-of-return  heuristic. 

Data  for  FY04  yields  491  MOSZs  that  are  eligible  for  an  SRB  multiplier.  With 
up  to  11  multiplier  values  for  some  MOSZs,  this  results  in  a  model  with  5,401  0-1 
variables  and  491  constraints.  A  solution  is  reached  in  1.4  seconds  on  a  1.58  GHz 
personal  computer.  The  solution  is  within  0.0018%  of  optimality  and  consumes  99.99% 
of  the  budget. 
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I.  INTRODUCTION 


Each  year  the  United  States  Marine  Corps  (USMC)  must  set  the  Selective 
Reenlistment  Bonus  (SRB)  for  each  Military  Occupational  Specialty  (MOS)  and  Service 
Time  Zone  (Zone)  (USMC  1990).  The  goal  is  to  selectively  encourage  Marines  to 
reenlist  and  enable  the  USMC  to  reach  reenlistment  targets  for  the  given  fiscal  year.  The 
purpose  of  this  thesis  is  to  give  manpower  planners  at  Manpower  &  Reserve  Affairs 
(M&RA),  Headquarters  Marine  Corps,  a  tool  that  will  identify  the  optimal  “SRB 
multiplier,”  for  each  MOS  and  Zone  combination  (MOSZ).  Theses  multipliers  set  the 
monetary  values  of  the  bonuses. 

A.  ENLISTED  MARINES,  MOS  AND  ZONE 

The  Marine  Corps  uses  two  models,  the  First  Term  Alignment  Plan  (FTAP)  and 
the  Subsequent  Term  Alignment  Plan  (STAP),  to  determine  the  number  of  reenlistments 
that  are  necessary  to  meet  the  future  force  structure  as  depicted  by  the  Grade  Adjusted 
Recapitulation  (GAR)  (Manpower  Plans  and  Policy  Branch-Integration  and  Analysis 
2007).  The  GAR  is  a  planning  tool  that  specifies  target  enlistment  numbers  for  all  MOSs 
and  ranks.  For  planning  purposes,  a  Marine  is  classified  by  pay  grade  and  time  in 
service.  Pay  grade  is  associated  with  a  Marine’s  level  of  responsibility  while  time  in 
service  correlates  to  experience.  A  Private  (E-l)  is  the  lowest  enlisted  pay  grade,  is  the 
least  experienced  Marine  with  less  than  one  year  in  service,  and  is  generally  still  in  initial 
training.  A  Sergeant  Major  (E-9)  is  the  highest  enlisted  pay  grade  and  is  the  most 
experienced  enlisted  Marine,  having  15  or  more  years  in  service. 

For  reenlistments  purposes,  manpower  planners  are  concerned  with  Marines  who 
are  in  or  about  to  enter  the  career  force,  meaning  they  have  21  or  more  months  time  in 
service.  These  Marines  are  placed  into  one  of  three  Service  Time  Zones  according  to 
time  in  service  as  listed  in  Table  1.  The  FTAP  model  covers  Marines  in  Zone  A  while 
the  STAP  model  covers  Marines  in  Zones  B  and  C. 
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Table  1.  Service  Time  Zones  designated  by  time  in  service. 


ZONE 

FROM 

TO 

A 

21  Months 

6  Years 

B 

6  Years 

10  Years 

C 

10  Years 

14  Years 

FTAP  is  the  retention  model  used  to  help  determine  the  number  of  reenlistments 
needed  from  Zone -A  Marines,  according  to  MOS,  to  meet  the  career-force  requirements 
set  by  the  GAR  (USMC  2004).  A  Zone-A  Marine  is  in  his  initial  enlistment  contract, 
including  extensions.  The  purpose  of  the  FTAP  is  to  fill  vacancies  in  the  career  force  that 
result  from  separation  or  loss  for  any  reason  while  preventing  promotion  stagnation  and 
ensuring  opportunities  for  advancement. 

STAP  is  similar  to  the  FTAP,  but  focuses  on  subsequent  reenlistments  of  Marines 
who  constitute  the  career  force  (Manpower  Plans  and  Policy  Branch-Integration  and 
Analysis,  2007).  The  primary  purpose  of  the  STAP  is  to  “move”  career-force  inventory 
levels  towards  required  levels.  STAP  recognizes  that  any  career-force  imbalances  cannot 
be  “fixed”  in  one  year,  and  tries  to  make  the  following  year  better  in  relation  to  the 
imbalance  in  the  previous  year. 

FTAP  is  important  because  it  covers  the  most  flexible  opportunity  the  Marine 
Corps  has  to  shape  the  career  force,  with  respect  to  the  number  of  Marines  in  each  MOS 
and  pay  grade.  Only  during  a  Marine’s  first  reenlistment  opportunity  can  the  Marine 
Corps  deny  reenlistment  (without  monetary  cost)  because  his  particular  MOS  is  full:  A 
reenlistment  is  allowed  only  when  there  is  a  specific  vacancy  in  the  career  force  that 
needs  to  be  filled. 

Provided  that  a  first-term  Marine  is  qualified,  any  Marine  who  is  denied 
reenlistment  in  his  original  MOS,  because  no  vacancies  exist,  is  allowed  to  reenlist  in 
another  MOS  that  needs  additional  manpower  (USMC  1993).  This  process  is  called  a 
Lateral  Move  or  “Lat  Move.”  Reenlistments  for  subsequent-tenn  (Zone-B  and  Zone-C) 
Marines  are  not  based  on  vacancies  but  instead  are  based  on  targets.  Furthermore,  the 
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MOS  will  not  close  to  reenlistments  even  if  the  target  is  reached.  Although  allowed  in 
specific  cases,  Lat  Moves  are  not  used  with  subsequent-tenn  Marines  as  they  are  with 
first-term  Marines. 

The  ability  to  deny  first-term  reenlistment  in  conjunction  with  the  use  of  Lat 
Moves  allows  manpower  planners  to  control  the  number  of  Marines  that  advance  into  the 
career  force  from  the  first-term  population,  in  specific  MOSs,  and  to  maintain  the 
personnel  goals  set  forth  by  the  GAR.  During  subsequent  reenlistment  opportunities  with 
a  career-force  Marine,  the  Marine  cannot  be  denied  reenlistment  provided  the  Marine  has 
no  performance  or  criminal-conduct  problems  and  has  been  promoted  on  schedule.  If,  in 
order  to  control  the  MOS  population,  manpower  planners  need  to  separate  a  subsequent- 
tenn  Marine  who  has  no  difficulties,  the  affected  Marine  is  entitled  to  a  substantial 
involuntary-separation  allowance . 

Historically,  most  MOSZs  do  not  need  a  bonus  to  encourage  enough  Marines  to 
reenlist  to  meet  assigned  targets.  Approximately  half  of  the  bonus  eligible  MOSZs  will 
not  receive  a  bonus  because  the  anticipated  manpower  shortfall  is  not  severe  enough  to 
warrant  a  bonus.  For  the  remaining  MOSZs,  the  manpower  shortfall  is  critical  and  the 
only  option  that  can  help  meet  the  remaining  reenlistment  requirements  is  the  SRB. 

B.  HOW  THE  SRB  IS  USED 

When  too  few  Marines  want  to  reenlist  or  Lat  Move  into  a  critical  MOSZ,  the 
Marine  Corps  uses  the  SRB  to  increase  the  number  of  reenlistments  there.  Typically,  a 
critical  MOS  falls  into  one  of  the  following  categories  (Morgan  2006): 

•  Technical  skills  requiring  large  training  and/or  replacement  costs, 

•  Skills  in  high  demand  in  the  civilian  sector, 

•  Skills  for  which  recruiting  is  difficult, 

•  Skills  that  are  crucial  to  combat  readiness,  and 

•  Skills  in  high  demand  but  with  a  small  reenlistment  population. 
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The  Marine  Corps  separates  Marines  into  one  of  three  Service  Time  Zones  by 
their  time  in  service,  as  defined  in  the  previous  section.  An  SRB  can  be  awarded  to  a 
Marine  only  once  in  each  Zone,  but  a  Marine  may  be  eligible  to  receive  an  SRB  in  each 
Zone  through  his  career. 

The  monetary  amount  awarded  for  an  SRB  is  computed  as  a  function  of  the 
Marine’s  pay,  additional  years  of  service  that  reenlistment  obligates,  and  the  “strength”  of 
the  bonus  referred  to  as  an  SRB  multiplier.  More  precisely,  the  SRB  is  defined  as 

The  Minimum  of: 

Monthly  Base  Pay  x  Years  of  Additional  Obligated  Service 
x  SRB  Multiplier  for  MOSZ 
and 

Maximum  Allowable  Bonus  for  the  Fiscal  Year. 

The  Marine  Corps  uses  SRB  multipliers  ranging  from  0  to  5  in  increments  of  0.5. 
In  an  MOS  that  does  not  need  an  incentive;  the  multiplier  is  0  and  no  bonus  is  offered.  In 
an  MOS  with  an  expected  shortfall,  the  multiplier  will  typically  be  non-zero. 

For  Zone  A,  manpower  planners  have  the  power  to  control  the  amount  of  total 
bonus  money  that  is  paid  by  closing  an  MOS  to  reenlistment.  Manpower  planners  cannot, 
however,  close  off  bonuses  or  reenlistments  to  a  Zone-B  or  C  MOS,  once  a  bonus  is 
offered  there.  This  means  that  if  all  of  the  eligible  Marines  in  an  MOS  reenlist,  all  will 
receive  the  bonus  money.  These  reasons  combined  with  the  limited  yearly  budget  for  the 
SRB  Program  make  picking  the  appropriate  multiplier  values  crucial. 

C.  CURRENT  SRB  MODELS 

Currently,  the  SRB  Program  uses  two  manual  heuristics  to  specify  SRB  multiplier 
values,  one  heuristic  for  FTAP  (Zone  A)  and  the  other  for  STAP  (Zones  B  and  C).  The 
total  budget  is  first  divided  between  FTAP  and  STAP,  and  then  the  heuristics  are  run 
separately.  Both  heuristics  are  implemented  using  Excel  spreadsheets  that  first  set  the 
multiplier  in  each  MOSZ  at  the  minimum  level  that  is  estimated  to  lead  to  the  meeting  of 
its  reenlistment  requirement.  In  order  to  detennine  which  multiplier  to  assign,  both 
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models  use  expected  reenlistment  rates  provided  by  the  Center  for  Naval  Analysis  (CNA) 
(North  1995).  Once  the  initial  settings  are  complete,  the  expected  costs  for  the  selected 
multipliers  are  summed  over  all  MOSZs  to  estimate  a  total  cost  of  the  program.  If  this 
estimate  exceeds  the  SRB  budget,  individual  SRB  multipliers  for  less  critical  MOSs  are 
manually  reduced  until  it  does  not. 

Manpower  planners  combine  the  output  of  FTAP  and  STAP  into  a  list  of 
prospective  SRB  multipliers,  and  carry  out  an  iterative,  negotiating  process  with  other 
departments  at  Headquarters  Marine  Corps,  to  agree  upon  the  final  SRB  multiplier 
values.  The  other  departments  include  the  occupational-field  sponsors,  assignment 
monitors,  command-level  representatives,  and  career-retention  specialists,  all  of  whom 
are  the  end-users  of  the  SRB  Program. 

D.  SUMMARY  OF  PAST  WORK  ON  SRB  OPTIMIZATION  AND  OTHER 

RELATED  RESEARCH 

A  1986  Naval  Postgraduate  School  thesis  by  U.S.  Marine  Corps  Captain  Dean  D. 
DeWolfe,  entitled  Determination  of  Selective  Reenlistment  Bonus  Multipliers  in  the 
United  States  Marine  Corps,  is  the  basis  for  this  thesis.  First,  DeWolfe  sets  up  and 
defines  the  problem  of  selecting  an  SRB  multiplier  for  each  MOSZ  as  a  nonlinear 
knapsack  problem.  The  objective  is  to  select  multipliers  that  minimize  a  total  weighted 
deviation  from  the  number  of  estimated  reenlistments  versus  the  reenlistment  targets, 
subject  to  a  budget  constraint;  all  coefficients  are  expressed  as  expected  values.  Next, 
DeWolfe  reformulates  the  problem  as  a  special  generalized  assignment  problem  shows 
how  to  maximize  a  Lagrangian  lower  bound  on  the  optimal  objective  value  and,  in  the 
process,  identify  a  good  feasible  solution.  Finally,  due  to  the  nature  of  the  Lagrangian 
relaxation,  it  is  possible  the  solution  does  not  use  the  entire  budget.  In  that  case  a 
secondary  heuristic  tries  to  increase  multipliers  to  consume  it  all.  The  heuristic  uses  a 
marginal  rate  of  return  to  detennine  which  MOS  multiplier  value  can  be  raised  to  the  next 
level  to  achieve  the  largest  decrease  in  the  objective  value  per  unit  of  (additional)  budget 
consumed. 
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Although  M&RA  liked  how  the  model  solved  the  SRB  multiplier  problem,  there 
were  problems  that  made  it  unusable.  First,  the  model  was  wrihen  in  FORTRAN  77 
which  is  not  familiar  to  personnel  at  M&RA.  Second,  there  have  been  many  changes  to 
the  SRB  program  rules  since  1986  that  needed  to  be  updated  (SRB  multiplier  increments, 
budgets,  lump  sums,  maximum  allowed  bonuses),  and  many  of  these  rules  were  “hard¬ 
wired”  into  the  code.  Finally,  and  most  importantly,  M&RA  wanted  a  model  that  could 
run  in  a  familiar  computational  environment  like  Excel’s  VBA  (Visual  Basic  for 
Applications).  Such  an  environment  would  allow  them  to  internally  use  and  maintain  the 
model. 

Other  related  research  includes  two  NPS  theses  by  U.S.  Marine  Corps  Captain 
Jonathan  D.  Raymond  (2006)  entitled  Determining  the  Number  of  Reenlistments 
Necessary  to  Satisfy  Future  Force  Requirements,  and  by  U.S.  Marine  Corps  Major  Dean 
G.  Conatser  (2006)  entitled  Forecasting  U.S.  Marine  Corps  Reenlistments  by  Military 
Occupational  Specialty  and  Grade.  Raymond  created  a  model  to  detennine  the  number 
of  reenlistments  that  are  necessary  for  each  MOS  and  grade  by  applying  transition  rates 
to  the  current  inventory  of  Marines  who  are  not  eligible  for  reenlistment  and  subtracting 
that  amount  from  the  desired  force  structure  in  the  GAR.  Conatser  created  a  model  to 
predict  the  number  of  reenlistments  for  both  FTAP  and  STAP  Marines  by  MOS  and 
grade.  These  two  models,  developed  at  the  request  of  M&RA,  are  the  beginning  of  the 
development  of  the  Career  Force  Retention  Model  (CFRM)  which  is  being  explored  to 
replace  the  FTAP  and  STAP  models.  The  CFRM,  in  conjunction  with  this  thesis,  will 
help  manpower  planners  estimate  the  required  number  of  new  recruits,  reenlistments,  and 
the  most  effective  way  to  set  SRB  multipliers  to  influence  the  reenlistments. 

This  thesis  extends  and  updates  DeWolfe’s  model  and  solution  algorithm  for  use 
by  M&RA.  This  is  accomplished  by  revising  budget  constraints,  adding  a  Lat  Move 
functionality,  and  making  all  model  parameters  accessible  to  the  user  (i.e.,  no  hard-wiring 
of  parameter  values  into  the  code).  The  ultimate  goal  of  this  thesis  is  to  create  a  model 
that  is  user-friendly  and  “user-updatable.” 

In  a  paper  by  DeWolfe,  Stevens,  and  Wood  in  1993  entitled  Setting  Military 

Reenlistment  Bonus,  the  authors  expand  upon  the  model  that  was  developed  in 
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DeWolfe’s  1986  thesis  for  use  by  the  United  States  Army.  First  a  more  complicated 
formulation  was  created  by  adding  an  additional  constraint  that  limited  the  number  of 
“high-level”  bonuses  to  10%  of  the  total  bonuses  offered.  Second,  the  formulation  found 
multipliers  for  each  MOS  and  Zone  based  on  expected  reenlistment  rates  for  3,  4,  5,  and 
6  years  of  additional  obligated  service.  Since  the  10%  constraint  is  no  longer  current  and 
the  Marine  Corps  only  offers  four-year  reenlistments,  the  more  complicated  formulation 
is  not  necessary  in  this  thesis. 

E.  THESIS  OUTLINE 

In  Chapter  II,  the  SRB  multiplier  problem  is  formulated  as  a  generalized 
assignment  problem.  The  data  and  variables  are  described  along  with  the  methodology 
used  to  compute  model  coefficients  like  the  training  costs  in  the  object  function.  This 
model  is  then  implemented  and  solved  in  the  General  Algebraic  Modeling  System 
(“GAMS”;  see  Rosenthal  2007),  with  an  Excel  interface,  for  testing  purposes.  Although 
the  GAMS  implementation  and  solution  is  a  perfectly  reasonable  one,  it  is  unsatisfactory 
for  the  final  product  of  this  thesis,  because  it  relies  on  an  unfamiliar  computational 
environment  for  manpower  planners.  Instead,  this  solution  is  used  in  Chapter  III  to  check 
the  accuracy  of  a  solution  method  based  on  Lagrangian  relaxation,  whose  implementation 
and  testing  in  VBA  and  Excel  is  also  described  in  Chapter  III.  Chapter  IV  presents 
conclusions  and  recommendations  for  further  work.  Appendices  A  and  B  contain  the 
source  code  for  the  GAMS  and  Excel  VBA  models,  respectively.  Appendix  C  contains 
sample  inputs  and  outputs  from  the  Excel  VBA  model. 
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II.  GENERALIZED  ASSIGNMENT  MODEL 


This  chapter  models  the  optimization  of  SRB  multipliers  across  all  MOSZs  as  a 
generalized  assignment  problem  (GAP).  First,  the  model  terms  and  coefficients  are 
defined.  Second,  the  differences  from  DeWolfe’s  model  and  the  current  model  are 
explained.  Finally,  the  estimation  methods  used  for  some  of  the  coefficients  in  the  model 
are  defined. 

A.  PROBLEM  FORMULATION 

There  is  a  total  of  m  MOSZs  that  encompass  all  MOSs  and  Zones.  Each  MOSZ 
will  take  on  one  of  n  multiplier  “levels,”  with  level  j  =  1  corresponding  to  a  multiplier 
value  of  Vj=0.  Current  SRB  policy  uses  multipliers  in  increments  of  0.5,  up  to  a 
maximum  of  5.0,  for  a  total  of  1 1  multiplier  levels.  Thus,  each  increase  in  level  j  results 
in  an  increase  of  0.5  in  v, .  For  example,  level  j  =  2  corresponds  to  v2  =  0.5  ,  j  =  3  to 

v3  =1.0,  and  the  maximum  level  of  j  =  11  corresponds  to  vu  =5.  The  model  is  not 
limited  to  a  predefined  number  of  multiplier  levels,  as  long  as  v  .+1  -  v .  is  constant  for  all 
levels  j . 

Currently,  there  is  no  set  maximum  multiplier  value  assigned  to  any  single  MOSZ 
or  group  of  MOSZs  that  is  different  from  the  overall  maximum  multiplier  mentioned 
above.  However,  the  multiplier  cannot  exceed  a  value  that  would  result  in  a  bonus  that 
exceeds  the  maximum  allowable  bonus.  The  maximum  allowable  multiplier  is  easy  to 
compute  and  is  represented  through  ni ,  the  maximum  multiplier  level  for  MOSZ  i . 

Each  MOSZ  has  an  associated  cost  (budget  consumption)  and  a  penalty 
(objective-function  value)  that  depends  on  the  expected  number  of  reenlistments  for  the 
chosen  multiplier.  The  overall  objective  is  to  assign  each  MOSZ  a  multiplier  that 
minimizes  the  total  penalty  and  keeps  the  total  cost  under  budget.  The  following 
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formulation  of  the  SRB  multiplier  problem  as  a  GAP  is  similar  to  DeWolfe’s  but  some 
variables  and  parameters  have  been  changed  to  account  for  Lat  Moves  and  the  use  of 
non-integer  multipliers. 

1.  Indices 

i  =  MOSZ 

j  =  1 allowable  SRB  multiplier  levels  for  MOS  i 

2.  Parameters  [units] 

d ..  reenlistment-target  deviation  penalty  for  MOSZ  i  with 

multiplier  level  j  [arbitrary  penalty  units] 

c..  budget  outlay  incurred  by  setting  multiplier  level  to  j  for  MOSZ  i  [dollars] 

B  budget  [dollars] 

3.  Decision  Variables 

Xj  1  if  multiplier  level  j  is  selected  for  MOSZ  i,  0  otherwise 

4.  Formulation  for  GAP 

m  nj 

z*  =min  (PI) 

i= 1  7=1 

m  Jij 

s.t  (Cl) 

1=1  7=1 

ZV=1  VI  (C2) 

7=1 

x,e{0,l  }V/,y  (C3) 

The  objective  function  (PI)  sums  penalties  across  all  MOSZs.  The  penalty  is  0 

for  a  reenlistment  target  that  is  met  exactly,  so  the  best  possible  objective-function  value 
is  0.  The  first  constraint  (Cl)  limits  total  bonus  expenditures  to  the  available  budget  (at 
least  in  terms  of  expected  values).  Note  that  cn  =  0  for  all  MOSZs  i ,  that  is,  no  cost  is 
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incurred  if  no  bonus  is  offered.  The  second  constraint  (C2)  ensures  that  each  MOSZ 
receives  exactly  one  multiplier.  This  model  follows  DeWolfe’s  exactly,  except  that 
coefficients  d  and  c  are  computed  differently,  as  described  below: 


T 


W, 


d=MTC  x 
v  A 


max 


{0, (/J  - LMiLMFi )  - RRRjE, ! 

+  Q  x  max  {o,  RRRijEi  -  (/J  -  LMiLMFi  )}2 


(Eqn  2.1) 


O  RRR  F  FvT  (Eqn  2.2) 

where  coefficients  are  defined  by 

Dj  number  of  Marines  needed  to  reenlist  in  MOSZ  i  [Marines] 

Wj  exogenous  weighting  factor  for  MOSZ  i  [unitless] 

P  average  monthly  pay  for  Marines  in  MOSZ  i  [dollars] 
v .  numerical  multiplier  value  for  multiplier  level  j 

RRR,,  expected  reenlistment  fraction  for  Marines  in  MOSZ  i  and  multiplier  level  j  [fraction] 
B  budget  [dollars] 

Q  relative  penalty  weight  of  overages  to  underages  [fraction] 

Y  average  years  of  reenlistment  [years] 

MTC  maximum  training  cost  for  all  MOSZs  [dollars] 

LMf  number  of  MOSZ  i  school  seats  open  for  Lateral  Moves  [school  seats] 

LMFj  fraction  of  available  school  seats  to  use  for  Lateral  Moves  into  MOSZ  i  [fraction] 

7]  training  cost  of  MOSZ  i  [dollars] 

At  total  number  of  Marines  in  MOSZ  i  [Marines] 

E.  number  of  Marines  eligible  for  reenlistment  in  MOSZ  i  [Marines] 


In  Eqn  2.1,  a  Lat-Move  functionality  is  added  by  the  term  LM iLMFj .  This  term 
can  reduce  the  number  of  reenlistments  that  is  needed  from  the  eligible  group  of  Marines 
in  each  MOSZ.  Each  MOSZ  has  a  predetermined  number  of  MOS  school  seats,  LMt , 

that  are  open  for  Lat  Moves.  Depending  on  the  MOSZ,  the  number  of  school  seats  can 
range  from  zero  to  the  total  number  of  Marines  that  are  needed  to  reenlist.  In  some 
MOSZs,  manpower  planners  will  want  to  use  all  available  school  seats  and  set 
LMFj  =1.0.  In  others,  planners  will  want  more  control  over  the  number  of  new  Marines 
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that  LatMove  into  the  MOSZ  and  set  LMF:  <  1 .  Manpower  planners  can  control  the 
number  of  LatMoves  by  assigning  LMFj ,  a  fraction  of  the  school  seats  to  use  for  each 

MOSZ.  Historically  there  are  always  enough  Lat-Move  eligible  Marines  to  fill  school 
seats  in  a  given  FY,  but  manpower  planners  will  need  to  verify  that  this  assumption  is 
still  valid. 

There  are  two  changes  to  Eqn  2.2  over  DeWolfe’s  formulation.  First,  the  total 
bonus  amount  for  each  MOSZ  is  used  in  the  budget  constraints  to  account  for  the  use  of 
lump-sum  payments  instead  of  the  anniversary  payments  in  use  in  1986.  Second,  a  new 
term  v,  is  introduced  to  account  for  the  use  of  non-integer  SRB  multipliers.  It  is  a 

parameter  for  the  SRB  multiplier  value  that  is  used  in  calculating  the  expected  cost  of 
using  an  arbitrary  multiplier  v . .  With  this  new  parameter,  manpower  planners  will  not 

be  limited  to  the  integer-only  multipliers  of  1986  or  the  multipliers  increments  of  0.5  in 
use  today. 

B.  TRAINING-COST  DATA 

A  key  term  in  the  weighting  function  for  the  objective  of  the  SRB  model  is  the 
training-cost  parameter,  which  weights  the  deviation  penalty  linearly  with  increasing 
training  cost:  If  a  shortfall  of  X  Marines  must  occur  in  some  MOSZ,  it  is  better  for  that 
to  occur  in  an  MOS  where  the  “replacement  cost,”  based  on  training  costs,  is  lower. 
Unfortunately,  the  Marine  Corps  no  longer  tracks  training-school  costs,  because  multiple 
sources  of  funding  makes  this  accounting  task  difficult.  Currently,  the  Marine  Corps  can 
only  supply  the  school  costs  for  non-Marine  Corps  schools  because  the  Marine  Corps 
pays  a  known  amount  for  each  Marine  to  attend.  In  order  to  overcome  the  lack  of  actual 
cost  data,  a  new  training-cost  estimate  is  developed  here,  based  on  Zone,  length  of  school 
or  schools  attended,  graduation  rate,  and  an  estimate  of  the  student’s  daily  pay.  All 
MOSZ  training-cost  parameter  estimates  use  a  variation  of  the  following  equation: 

Training-cost  estimate  = 

(Training  Days  x  Daily  Pay)  /  Graduation  Rate  (Eqn  2.3) 
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1.  Training-Cost  Estimate  for  Zone  A 


Zone-A  training  costs  are  broken  into  two  categories:  Initial  training  MOSs  and 
Lat-Move-only  MOSs.  Both  categories  use  Eqn  2.3  to  estimate  corresponding  training 
costs  but  use  different  daily  pay  amounts  to  reflect  the  average  pay  received  by  a  Marine 
while  in  training.  This  thesis  applies,  for  initial  training  MOSs,  the  daily  pay  for  an  E-l 
with  less  than  two  years  in  service;  for  Lat-Move-only  MOSs,  it  uses  the  daily  pay  for  an 
E-5  with  more  than  four  years  in  service. 

2.  Training-Cost  Estimate  for  Zones  B  and  C 

Zone-B  and  -C  training  costs  are  also  based  on  Eqn  2.3  but  are  more  complicated 
because  of  the  various  tracks  that  a  Marines’  career  can  take  in  his  MOS,  an  MOS  that 
may,  in  fact,  change.  First,  the  Marine  can  be  promoted  and  keep  his  original  MOS. 
Under  those  circumstances,  training  cost  is  based  on  the  MOS’s  Zone-A  training  cost. 
Second,  a  Marine  can  be  promoted  and  receive  a  new  MOS  to  distinguish  new 
responsibilities.  The  new  MOS  can  be  supplied  by  a  single  MOS  or  by  multiple  “feeder 
MOSs.”  In  the  fonner  case,  the  training-cost  estimate  is  also  based  on  the  Zone-A 
training  cost  of  the  old  MOS;  otherwise  the  training  cost  is  the  average  total  training  costs 
of  the  Zone-A  feeder  MOSs.  In  all  of  the  above  cases,  it  is  possible  that  the  Marine  will 
be  required  to  attend  another  MOS  training  school.  If  applicable,  the  new  training  cost 
for  the  school  is  estimated  through  Eqn  2.3  and  is  added  to  the  prior-training  cost 
estimate  that  is  determined  by  the  MOS  track. 

C.  GAMS  MODEL 

Using  FY2004  data,  the  GAP  is  implemented  and  solved  in  GAMS  using  the  XA 
solver  (GAMS  2007).  The  FY2004  data  set  contains  all  809  USMC  MOSZs,  of  which 
491  MOSZs  need  a  multiplier  greater  than  zero  to  reach  their  reenlistment  goals.  The 
final  objective  value  is  57.8525  is  obtained  in  0.02  seconds  of  solver  time,  and  is  proven 
to  be  within  0.0006%  of  optimality.  Of  the  491  MOSZs,  385  receive  a  multiplier  greater 
than  zero.  Table  2  shows  the  distribution  of  the  eligible  MOSZs  and  their  selected 
multipliers. 


13 


Table  2.  Distribution  of  SRB  Multipliers  in  GAMS  Solution 


Multiplier 

0.0 

0.5 

1.0 

1.5 

2.0 

2.5 

3.0 

3.5 

4.0 

4.5 

5.0 

MOSZ 

106 

39 

99 

41 

33 

95 

74 

0 

1 

3 

0 

The  eligible  MOSZs  have  a  reenlistment  goal  of  8,887  Marines,  and  with  the 
chosen  multipliers  5,620.75  Marines  are  expected  to  reenlist.  The  SRB  budget  for 
FY2004  is  $54  million  and  only  $474  goes  unspent  in  the  GAMS  solution. 
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III.  SOLVING  THE  GENERALIZED  ASSIGNMENT  MODEL 
WITH  L AGRAN GIAN  RELAXATION 


Although  the  updated  GAP  formulation  of  the  SRB  multiplier  problem  solves 
quickly  in  GAMS,  manpower  planners  cannot  readily  use  this  implementation  due  to  lack 
of  access  to  GAMS.  In  order  to  supply  manpower  planners  with  a  model  they  can  use, 
the  GAP  model  is  implemented  and  solved  here  using  the  Lagrangian-relaxation 
technique  described  by  DeWolfe,  but  using  a  current  programming  language  that  can  run 
in  a  computational  environment  that  manpower  planners  have  access  to.  The  following 
section  borrows  heavily  from  DeWolfe  (1986),  and  the  reader  should  see  that  thesis  for 
more  details. 

A.  L  AGRAN  GIAN  RELAXATION 

The  budget  constraint  (Cl)  in  the  GAP  can  be  viewed  as  a  “complicating 
constraint.”  That  is,  the  model  solves  easily  if  (Cl)  is  relaxed.  Using  Lagrangian 
relaxation,  the  budget  constraint  can  be  moved  to  the  objective  function  using  a 
Lagrangian  multiplier  that  acts  as  a  penalty  for  constraint  violation,  or  as  an  incentive  to 
spend.  The  relaxed  formulation  is: 

For  A  >  0  , 

m  11 ,  in  " , 

LR(T) :  z(X)  =  min^  £  dyxy  +  AC£^ c..x..  -  B)  (P3) 

i=\  j= 1  i=l  7=1 

ni 

s.t.  YJxj=\  Vi 

)= i 

Xu  e  {0,1}  Vi,  j 

The  inner  portion  of  P3  can  be  simplified  and  rewritten  as: 

m  n  j 

minXZ Ky  +  Acij )x(/  ~AB  (Eqn  3.1) 

i=l  7=1 
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Thus,  for  fixed  X ,  z(X)  is  computed  by  solving  a  simple  selection  problem:  For 
each  MOSZ  i ,  choose  x..  to  be  1  for  the  largest  value  of  djj  +  /lev  ,  and  set  other  x..  to  0. 

It  is  well  known,  then,  that  z(X)  <  z * ,  i.e.,  z(X)  provides  a  lower  bound  on  z* . 

The  best  lower  bound  is  found  by  maximizing  z(/t)with  respect  to  X  to  obtain 
X* .  This  can  be  done  by  bracketing  the  optimal  X  in  some  range  \LnLu\,  and  then 
narrowing  the  bracket  through  bisection  search.  L,  =  0  is  always  valid.  For  X 
sufficiently  large,  the  solution  to  LR(/t)  will  spend  as  few  SRB  budget  dollars  as 
possible.  This  will  occur  when  setting  x, .  =1  for  all  i  is  the  optimal  action,  and  this  will 
occur  when 

da  +  Xcn  <  dj  +  Xcy  Vi,  j>  2  (Eqn  3 .2) 

Since  ca  =  0  ,  this  leads  to 

da^dy+Xcy  V/,/>2. 

Thus 

Lu  =  max  {{dn  -dv)/  ctj }  V/,  j  >  2 

is  a  valid  value  for  the  right-hand  side  of  the  bracket. 

The  bisection  search  is  based  on  the  fact  that  if  the  budget  constraint  is  violated 
for  a  given  X ,  then  X*  >  X  ;  otherwise,  X*  <  X  . 

The  above  methodology  is  implemented  just  as  in  DeWolfe  (1986),  but  using 
VBA,  and  the  best  Lagrangian  lower  bound,  z(X*)  ,  is  established.  An  upper  bound  on 
the  optimal  solution  to  the  SRB  multiplier  problem  is  found  by  examining  the  feasible 
solution  sets  to  z(X)  that  are  encountered  in  the  course  of  maximizing  z(X) .  (When 

X  >  X* ,  the  solution  to  LR(Z)  is  guaranteed  to  be  feasible  with  respect  to  the  budget 
constraint.) 


(Eqn  3.3) 

(Eqn  3.4) 
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The  best  solution  identified  while  optimizing  z(A)  may  not  consume  the  entire 
budget.  If  unspent  budget  remains,  it  may  be  possible  to  improve  the  solution  using  the 
marginal-rate-of-retum  heuristic  described  by  DeWolfe.  The  marginal  rate  of  return  for 
MOSZ  i  is  defined  as: 

R()R,  =  (dj  -  di  j+l )  /(ciJ+l  -  cy  )  (Eqn  3.5) 

where  the  current  feasible  solution  has x.  =  1 .  The  numerator  in  Eqn  3.5  calculates  the 
amount  z(X)  can  be  improved  by  increasing  the  SRB  level  for  MOSZ  i  from  its  current 
level  of  j  to  j  +  1;  the  denominator  calculates  the  corresponding  additional  budget 
expenditure.  In  essence,  RORj  calculates  the  “bang  for  the  buck”  for  raising  the 

multiplier  level  for  each  MOSZ  individually.  The  heuristic  repeatedly  applies  these  rules 
until  no  additional  budget  can  be  feasibly  consumed: 

a.  Given  the  current  values  for  jc..  ,  compute  RORj  for  each  MOSZ  i. 

b.  Select  MOSZ  i  with  the  largest  value  of  RORi  such  that  (a)  setting 
Xy  =  0  and  xt  .+1  =  1  does  not  violate  the  budget  constraint,  and  (b)  j  + 1 
does  not  exceed  the  largest  allowable  SRB  level  for  the  MOSZ. 

c.  Set  Xy  =  0  and  xt  .+1  =  1 . 

B.  EXCEL  VBA  MODEL  SOLUTION  BY  LAGRANGIAN  RELAXATION 

Using  FY2004  data,  the  Lagrangian-based  solution  method  for  the  GAP,  as 
described  above,  is  implemented  in  EXCEL  using  VBA.  For  the  FY2004  data  described 
previously,  the  final  objective  value  of  57.8691  is  obtained  in  1.4  seconds,  and  is  proven 
to  be  within  0.0018%  of  optimality.  Of  the  491  MOSZs  eligible  for  a  bonus,  385  receive 
a  multiplier  greater  than  zero.  Table  2  shows  the  distribution  of  the  eligible  MOSZs  and 
their  selected  multipliers. 

Table  3.  Distribution  of  SRB  Multipliers  in  VBA  Solution 


Multiplier 

0.0 

0.5 

1.0 

1.5 

2.0 

2.5 

3.0 

3.5 

4.0 

4.5 

5.0 

MOSZ 

106 

39 

94 

42 

35 

97 

74 

0 

1 

3 

0 
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The  eligible  MOSZs  have  a  reenlistment  goal  of  8,887  Marines,  and  5,621.51 
Marines  are  expected  to  reenlist  with  the  chosen  multipliers.  Of  the  $54  million  budget, 
only  $574  goes  unspent. 

C.  COMPARISON  OF  DIRECT  AND  L AGRAN GIAN  SOLUTIONS 

The  VBA  and  GAMS  solution  methods  produce  very  similar  outputs,  with  the 
VBA  model  having  the  higher  objective-function  value  by  0.0166  units,  or  by  0.0287%. 
Of  the  491  MOSZs  that  are  eligible  for  an  SRB,  both  models  assign  non-zero  multipliers 
to  the  same  385  MOSZs.  Of  these  MOSZs,  the  VBA  solution  assigns  only  1 1  multiplier 
values  that  differ  from  the  GAMS  solution.  Of  these,  nine  multipliers  increase  by  0.5  in 
the  VBA  solution,  one  multiplier  decreases  by  0.5,  and  one  multiplier  increases  by  1.5. 

The  largest  difference  in  objective-function  values  results  from  the  lone 
multiplier  that  decreases  by  0.5.  The  decrease  raises  the  objective  value  by  0.1057  units 
or  0.18%  of  the  total  VBA  objective-function  value.  For  the  MOSZ  with  the  largest 
multiplier  difference  in  the  VBA  solution,  the  multiplier  is  raised  from  the  GAMS- 
solution  value  of  1.5  to  3.0.  Despite  this  large  difference,  the  impact  on  the  objective 
value  is  negligible:  The  objective  value  decreases  by  0.1 106  units  or  0.06%  of  the  total 
objective  value.  This  small  change  in  the  objective  value  results  from  the  fact  that  the 
relevant  MOSZ  contains  only  two  Marines,  and  so  even  the  worst  multiplier  assignment 
contributes  little  to  the  overall  objective-function  value  (penalty). 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


In  summary,  the  model  developed  in  this  thesis  provides  Marine  Corps  manpower 
planners  a  new  tool  to  help  manage  the  Selective  Reenlistment  Bonus  (SRB)  program 
(USMC  1990).  The  model’s  solution  specifies  “SRB  multipliers”  for  each  Military 
Occupational  Specialty  (MOS)  and  Service  Time  Zone  (Zone)  combination  that  define 
the  bonus  amount,  and  are  projected  to  encourage  reenlistments  that  will  (a)  meet  annual 
reenlistment  targets  as  best  possible,  and  (b)  not  exceed  the  available  SRB  budget.  The 
new  extended  SRB  optimization  model  and  solution  procedure  incorporates  key  aspects 
of  the  methodology  developed  in  DeWolfe  (1986),  but  makes  these  improvements: 

•  Gives  manpower  planners  the  ability  to  exploit  Lateral  Moves, 

•  Allows  a  more  flexible  set  of  multiplier  values  that  need  not  be  integers, 

•  Makes  budget  and  other  parameters  “user  updateable,”  and 

•  Implements  the  model  and  solution  procedure  in  a  familiar  computing 
environment  for  manpower  planners,  and  at  no  cost  (beyond  software  and 
hardware  that  these  planners  already  possess). 

The  new  model  does  have  its  limitations.  When  using  the  Lateral-Move 
functionality,  the  critical  assumption  about  the  availability  of  Marines  eligible  for  “Lat 
Moves”  needs  to  be  checked.  Also,  the  expected  number  of  reenlistments  for  a  given 
multiplier  is  based  on  a  forecasting  model  that  is  assumed  perfect. 

Two  areas  that  could  improve  the  effectiveness  of  the  model  warrant  further 

study: 

•  The  development  of  a  database  that  more  accurately  reflects  true  training 
costs,  rather  than  the  one  used  in  this  thesis,  which  estimates  those  costs 
based  on  training  days  and  average  pay, 
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•  The  expansion  of  the  model  to  better  forecast  the  expected  cost  of  a  given 
multiplier  by  breaking  down  the  Marines  in  a  MOSZ  according  to  pay 
grade. 
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APPENDIX  A:  GAMS  SOURCE  CODE 


This  appendix  contains  the  source  code  for  the  nonlinear  optimization  problem 
solved  in  GAMS  as  a  comparison  and  check  for  the  Excel  VBA  solution. 

$TITLE  SRB  Multiplier  Optimization 

* - DEFAULTS - 

SOFFUPPER  OFFSYMLIST  OFFSYMXREF 

OPTIONS 
LIMROW  =  0 
LIMCOL  =  0 
ITERLIM  =  500000 
RESLIM  =  100000 
SOLPRINT  =  OFF 
DECIMALS  =  2 
LP  =  XA 
RMIP  =  XA 
MIP  =  XA 
OPTCR  =  0.001 
OPTCA  =  0 


SONTEXT 

Original:  20  Jan  2007 

Author  :  Kent  A.  Robbins,  Jr. 

Class  :  Thesis 
Revised:  10  August  2007 

-Deleted  Zone  Loop  (Now  just  associated  with  MOS) 

-Added  read-in  data  function 

-Dropped  Bonus  limit  constraint  which  is  accounted  for 
in  the  data  that  is  received  from  the  Excel  VBA 
spreadsheet  that  writes  the  files. 

-Added  fixing  variables  for  Multipliers 
-Added  output  into  Excel  file 
-dropped  OPTCA  from  0.0  to  .001 
Description:  SRB  Multiplier  Optimization  Program  that  selects  the 
optimum  multiplier  to  minimize  the  squared  weighted 
deviation  of  the  differences  between  the  expected 
and  needed  number  of  reenlistments.  This  program 
is  run  in  shell  from  Excel’s  VBA. 
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$OFFTEXT 


* - INDICES  - 

SETS 
i  MOS 

$INCLUDE  C:\Documents  and  SettingsYKent  Robbins\My  Documents\Thesis\GAMS 
Thesis  Data  Files\MOSZoneList.dat 

m  Multiplier 

$INCLUDE  C:\Documents  and  Settings\Kent  Robbins\My  Documents\Thesis\GAMS 
Thesis  Data  Files\MultipleList.dat 


* - DATA  - 

Scalar 

Q  "relative  weight  of  overages  to  shortages  wrt  reenlistment  target" 

$INCLUDE  C:\Documents  and  SettingsYfCent  RobbinsYMy  DocumentsYThesisYGAMS 
Thesis  Data  FilesYQ.dat 
B  "Total  allowed  budget" 

$INCLUDE  CADocuments  and  SettingsYKent  RobbinsYMy  DocumentsYThesisYGAMS 

Thesis  Data  FilesYBudget.dat 

Yrs  "number  of  years  for  reenlistment" 

$INCLUDE  CADocuments  and  SettingsYKent  RobbinsYMy  DocumentsYThesisYGAMS 
Thesis  Data  FilesYYrs.dat 

maxTmCost  "maximum  training  cost  of  all  MOS's" 

$INCLUDE  CADocuments  and  SettingsYKent  RobbinsYMy  DocumentsYThesisYGAMS 
Thesis  Data  FilesYMaxTrnCost.dat 


Parameter  T(i)  "normalized  cost  of  traing  a  Marine  in  MOS.Zone(i)" 

$INCLUDE  CADocuments  and  SettingsYKent  RobbinsYMy  DocumentsYThesisYGAMS 
Thesis  Data  FilesYTrnCost.dat 


Parameter  A(i)  "number  of  Marines  in  MOS.Zone(i)" 

$INCLUDE  CADocuments  and  SettingsYKent  RobbinsYMy  DocumentsYThesisYGAMS 
Thesis  Data  FilesYMOSPop.dat 


Parameter  E(i)  "number  of  Marines  in  MOS.Zone(i)  eligible  for  reenlistment" 
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$INCLUDE  C:\Documents  and  Settings\Kent  Robbins\My  Documents\Thesis\GAMS 
Thesis  Data  Files\EASPop.dat 


Parameter  D(i)  "number  of  Marines  required  to  reenlist  in  MOS.Zone(i)" 

$INCLUDE  C:\Documents  and  Settings\Kent  Robbins\My  Documents\Thesis\GAMS 
Thesis  Data  Files\BoatSpace.dat 


Parameter  W(i)  "exogenous  weighting  factor  for  MOS.Zone(i)" 

$  INCLUDE  C:\Documents  and  Settings\Kent  Robbins\My  Documents\Thesis\GAMS 
Thesis  Data  Files\WgtFactor.dat 


Parameter  P(i)  "average  monthly  base  bay  of  a  Marine  in  MOS.Zone(i)" 

$  INCLUDE  C:\Documents  and  SettingsYfCent  Robbins\My  Documents\Thesis\GAMS 
Thesis  Data  Files\AvgPay.dat 


Parameter  MultiplierNumber(m)  "corresponding  multiplier  associated  with  set  (m)" 
$INCLEIDE  C:\Documents  and  Settings\Kent  RobbinsYMy  Documents\Thesis\GAMS 
Thesis  Data  Files\MultipleNumbers.dat 


Parameter  Preset(i)  "gives  the  value  of  the  preset  multiplier  for  MOS.Zone(i)" 
$INCLEIDE  C:\Documents  and  Settings\Kent  Robbins\My  Documents\Thesis\GAMS 
Thesis  Data  Files\PresetMult.dat 


TABLE 

RRR(i,m)  "reenlistment  response  rate  for  MOS.Zone(i)  with  Multiplier(m)" 
$ONDELIM 

$INCLUDE  C:\Documents  and  Settings\Kent  Robbins\My  Documents\Thesis\GAMS 

Thesis  Data  Files\ReenlistmentRates.dat 

$OFFDELIM 

Parameter  TEPDev(i,m)  "deviation  of  D(i)-RRR(i,m)*E(i)  or  0  (target-expected)"; 
TEPDev(i,m)=  max(  0,  D(i)-RRR(i,m)*E(i)  ); 

Parameter  ETPDev(i,m)  "deviation  of  RRR(i,m)*E(i)-D(i)  or  0  (expected-target)"; 
ETPDev(i,m)=  max(  0,  RRR(i,m)*E(i)-D(i)  ); 

Parameter  Wght(i)  "Weight  function  for  MOS.Zone(i)"; 

Loop((i),  IF(A(i)  >  0,  Wght(i)=  (T(i)/maxTrnCost)*W(i)/A(i); 
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ELSE  Wght(i)=0;);); 


*  - VARIABLE - 

BINARY  VARIABLE 

y(i,m)  Binary  var  is  1  if  Multiplier  m  is  used  for  MOS.Zone(i) 
VARIABLE 

Z  Total  deviation  for  selection  optimal  multipliers 

*  - EQUATIONS - 

EQUATIONS 

TOTALDEVIATION  total  cost  of  SRB  allocation 
PickOne(i)  only  allow  one  SRB  Multiplier 
TotalCost  total  cost  of  picking  mult  m  for  MOS.Zone(i) 


* - OBJECTIVE  FUNCTION - 

TOTALDEVIATION.. 

Z  =E=  Sum((i,m),(Wght(i)*(TEPDev(i,m)**2+Q*ETPDev(i,m)*Hs2))H5y(i,m)) 


* - CONSTRAINTS - 

PickOne(i).. 

Sum(m,  y(i,m))  =E=  1 ; 

TotalCost.. 

Sum((i,m),  RRR(i,m)*E(i)*y(i,m)*MultiplierNumber(m)*P(i)*Yrs)  =L=  B; 


*This  loop  will  fix  multipliers  if  the  MOS.Zone(i)  as  a  preset  values. 
LOOP((i)$(Preset(i)>-l), 

Loop((m),  if(MultiplierNumber(m)  =  Preset(i),  y.fx(i,m)  =  1;);); 

); 


MODEL  SRB  /ALL/ ; 

SOLVE  SRB  USING  MIP  MINIMIZING  Z  ; 

file  out  /CADocuments  and  Settings\Kent  RobbinsYMy  Documents\Thesis\GAMS  Thesis 
Data  Files'  SRBOptResults.csv/; 
put  out; 

*This  loop  writes  the  solution  to  a  file  that  will  be  read  in  by  Excel. 
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LOOP(i, 

LOOP(m$(y.l(i,m)=l), 
put  MultiplierNumber(m),’ 

); 

put/; 

); 


putclose  out; 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 
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APPENDIX  B:  VBA  SOURCE  CODE 


This  appendix  contains  the  source  code  for  the  VBA  modules  that  are  used  to  set 

up  the  SRB  optimization  problem  for  both  the  GAMS  and  VBA  models.  This  appendix 

also  contains  the  module  that  solves  the  Lagrangian  relaxation  to  the  Generalized 

Assignment  Model  solved  in  GAMS. 

'File:  Checks  Module  in  VBA 
'Date  Created:  12  December  2006 
'Last  Updated:  15  August  2007 

t 

'Checks  is  a  VBA  Module  that  contains  a  list  of  subroutines  that  will  perform  checks  on 
'the  data  set  to  verify  that  it  is  complete  in  order  to  run  the  optimization  software. 

t 

'@  Author  Kent  Robbins,  Jr. 

t 

Option  Explicit 
Public  INumNodes  As  Long 
Public  IMaxIndex  As  Long 
Public  sNodeNames()  As  String 
Public  llndex()  As  Long 
Public  INumErrors  As  Long 

'FillZoneABtSpaces()  is  a  subroutine  that  will  compute  the  correct  number  of 
'reenlistments  for  Zone  A  that  are  needed  when  BNA  school  seats  are  used. 

t 

'@param  -  None 

t 

Sub  LillZoneABtSpacesO 
Dim  rZoneABtSpacesAvail  As  Range 
Dim  rZoneABNASeatsAvail  As  Range 
Dim  rZoneABNASeatUsage  As  Range 
Dim  rZoneABtSpaceAdjust  As  Range 
Dim  iLength  As  Integer,  i  As  Integer 

Set  rZoneABtSpacesAvail  =  wsMOSZoneA.Range("ZoneABtSpaceStarf _ 
wsMOSZoneA.Range("ZoneABtSpaceStart").End(xlDown)) 

Set  rZoneABNASeatsAvail  =  wsMOSZoneA.Range("ZoneABNASeatStart",  _ 
wsMOSZoneA.Range("ZoneABNASeatStart").End(xlDown)) 

Set  rZoneABNASeatUsage  =  wsMOSZoneA.Range("ZoneABNASeatUsageStarf _ 
wsMOSZoneA.Range("ZoneABNASeatUsageStart").End(xlDown)) 
iLength  =  rZoneABtSpacesAvail.Rows.  Count 

Set  rZoneABtSpaceAdjust  =  wsMOSZoneA.Range("ZoneABSAdjStart",  _ 

wsMOSZoneA.Range("ZoneABSAdjStart").Offset(iLength,  0)) 
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For  i  =  1  To  iLength 

rZoneABtSpaceAdjust(i)  =  Application.WorksheetFunction.Max  _ 
((rZoneABtSpacesAvail(i)  -  Application.  _ 
WorksheetFunction.RoundUp(rZoneABNASeatsAvail(i)  *  _ 
rZoneABNASeatUsage(i),  0)),  0) 

Next  i 
End  Sub 

'MOSZoneCheckQ  is  a  subroutine  that  will  verify  that  the  MOS  is  listed  in  the  Training 
'Cost  Worksheet.  This  subroutine  gathers  the  data  and  calls  subroutine  CheckMOS  to 
'verify  that  MOS  is  listed.  If  there  are  MOSs  that  are  not  listed  a  window  will  pop  up 
'and  warn  the  user. 

t 

'@param  -  None 

t 

Sub  MOSZoneCheck() 

'Data  Sheet  Ranges 

Dim  rZoneAMOS  As  Range 

Dim  rZoneBMOS  As  Range 

Dim  rZoneCMOS  As  Range 

Dim  rZoneAMOSTrnVal  As  Range 

Dim  rZoneBMOSTmVal  As  Range 

Dim  rZoneCMOSTmVal  As  Range 

’Training  Sheet  Ranges 

Dim  rTrainMOSAList  As  Range 

Dim  rTrainMOSBCList  As  Range 

Dim  rTrainMOSACostList  As  Range 

Dim  rTrainMOSBCCostList  As  Range 

'Check  &  Fill  in  Zone  A  Ranges 

Set  rZoneAMOS  =  wsMOSZoneA.Range("ZoneAMOSStart",  wsMOSZoneA.Range 
(”ZoneAMOSStart").End(xlDown)) 

Set  rZoneAMOSTrnVal  =  wsMOSZoneA.Range("ZoneATrnCstStart",  _ 
wsMOSZoneA.Range("ZoneATrnCstStart").End(xlDown)) 

Set  rTrainMOSAList  =  wsTrainCostZoneA.Range("TrnMOSZoneAListStart",  _ 
wsTrainCostZoneA.Range("TrnMOSZoneAListStart")  _ 

.End(xlDown)) 

Set  rTrainMOSACostList  =  wsTrainCostZoneA.Range("TrnMOSZoneACostStart",  _ 
wsTrainCostZoneA.Range("TrnMOSZoneACostStart")  _ 
.End(xlDown)) 

Call  CheckMOS(rZoneAMOS,  rZoneAMOSTrnVal,  rTrainMOSAList,  _ 
rTrainMOSACostList) 

If  INumErrors  >  0  Then 

MsgBox  ("Error:  There  are  MOS(s)  in  Zone  A  not  listed  in  the  Training"  _ 

&  "  Cost  Data-Zone  A  worksheet."  &  vbCrLf  &  "  See  Data  -  "  _ 

&  "Zone  A  Workseet.") 
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End  If 

'Check  &  Fill  in  Zone  B  Ranges 

Set  rZoneBMOS  =  wsMOSZoneB.Range("ZoneBMOSStart",  wsMOSZoneB. Range  _ 
("ZoneBMOSStart").End(xlDown)) 

Set  rZoneBMOSTmVal  =  wsMOSZoneB.Range("ZoneBTmCstStart",  _ 
wsMOSZoneB. Range("ZoneBTmCstStart").End(xlDown)) 

Set  rTrainMOSBCList  =  wsTrainCostZoneBC.Range("TmMOSZoneBCListStart",  _ 
wsTrainCostZoneBC.Range("TrnMOSZoneBCListStart")  _ 
.End(xlDown)) 

Set  rTrainMOSBCCostList  =  wsTrainCostZoneBC. Range  _ 
("TrnMOSZoneBCCostStart",  _ 
wsTrainCostZoneBC.Range  _ 
("TrnMOSZoneBCCostStart").End(xlDown)) 

Call  CheckMOS(rZoneBMOS,  rZoneBMOSTmVal,  rTrainMOSBCList,  _ 
rTrainMOSBCCostList) 

If  INumErrors  >  0  Then 

MsgBox  ("Error:  There  are  MOS(s)  in  Zone  B  not  listed  in  the  Training"  _ 

&  "  Cost  Data-Zone  BC  worksheet."  &  vbCrLf  &  "  See  Data  -  "  _ 

&  "Zone  B  Workseet.") 

End  If 

'Check  &  Fill  in  Zone  C  Ranges 

Set  rZoneCMOS  =  wsMOSZoneC.Range("ZoneCMOSStart",  wsMOSZoneC. Range  _ 
("ZoneCMOSStart").End(xlDown)) 

Set  rZoneCMOSTmVal  =  wsMOSZoneC.Range("ZoneCTmCstStart",  _ 
wsMOSZoneC.  Range("ZoneCTmCstStart").End(xlDown)) 

Call  CheckMOS(rZoneCMOS,  rZoneCMOSTmVal,  rTrainMOSBCList,  _ 
rTrainMOSBCCostList) 

If  INumErrors  >  0  Then 

MsgBox  ("Error:  There  are  MOS(s)  in  Zone  C  not  listed  in  the  Training"  _ 

&  "  Cost  Data-Zone  BC  worksheet."  &  vbCrLf  &  "  See  Data  -  "  _ 

&  "Zone  C  Workseet.") 

End  If 
End  Sub 

'CheckMOSQ  is  a  subroutine  that  will  verify  that  the  MOS  is  listed  in  the 
'Training  Cost  Worksheet.  If  the  MOS  is  listed  the  training  cost  will  be  listed 
'in  the  training  cost  column,  if  not  an  error  will  be  written. 

f 

'@param  -  rMOSToCheck  is  the  range  of  MOSs  in  the  worksheet  that  need  to  be  verified 

-  rTrnValToFill  is  the  range  where  training  costs  will  be  filled  in  the 

'  worksheet 

-  rMOSList  is  the  range  of  MOSs  that  are  listed  in  the  Traing  Cost  Worksheet 

-  rTrnCost  is  the  range  of  training  costs  for  MOSs  listed  in  the  Traing  Cost 

'  Worksheet 
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Sub  CheckMOS(rMOSToCheck  As  Range,  rTrnValToFill  As  Range,  _ 
rMOSList  As  Range,  rTmCost  As  Range) 

Dim  lNumMOS  As  Long 
Dim  INumMOSToCheck  As  Long 
INumErrors  =  0 
INumNodes  =  0 

lNumMOS  =  rMOSList.Rows. Count 
INumMOSToCheck  =  rMOSToCheck.Rows. Count 
IMaxIndex  =  3  *  lNumMOS 
ReDim  llndex(l  To  IMaxIndex)  As  Long 
ReDim  sNodeNames(0  To  2  *  lNumMOS)  As  String 
Dim  i  As  Integer 
For  i  =  1  To  lNumMOS 
AddNode  rMOSList(i) 

Next  i 

rTmValToFill. Clear 

For  i  =  1  To  INumMOSToCheck 

If  DoesContain(rMOSToCheck(i))  Then 
rTmValToFill(i)  =  rTrnCost(i) 

Else 

rTmValToFill(i)  =  "Error,  MOS  not  LISTED" 

INumErrors  =  INumErrors  +  1 
End  If 
Next  i 

Call  F ormatT woDecimal(rT rnV alT oFill) 

End  Sub 

'ReenlistmentRateCheckByOccFldQ  is  a  subroutine  that  will  verify  that  the  OccField 
'is  listed  in  the  Reenlistment  Rate  Worksheet  for  the  correct  zone.  This 
'subroutine  will  call  ReenlistmentrateCheckMaster  and  make  it  verify  by  OccField. 

t 

'@param  -  None 

t 

Sub  ReenlistmentRateCheckByOccFld() 

Call  ReenlistmentRateCheckMaster(True) 

End  Sub 

ReenlistmentRateCheckByMOS()  is  a  subroutine  that  will  verify  that  the  MOS  is  listed 
'in  the  Reenlistment  Rate  Worksheet  for  the  correct  zone.  This  subroutine  will  call 
'ReenlistmentrateCheckMaster  and  make  it  verify  by  MOS. 

t 

'@param  -  None 

f 

Sub  ReenlistmentRateCheckByMOS() 

Call  ReenlistmentRateCheckMaster(False) 
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End  Sub 


ReenlistmentRateCheckMaster()  is  a  subroutine  that  will  verify  that  the  MOS/OccField 
'is  listed  in  the  Reenlistment  Rate  Worksheet  for  the  correct  zone.  This  subroutine 
'gathers  the  data  to  use  and  calls  the  subroutine  CheckRRRQ  to  verify  that  the 
'MOS/OccField  is  listed.  If  the  MOS/OccField  is  listed  the  reenlistment  rates  and  its 
'corresponding  multiplier  will  be  listed  in  the  reenlistment  rate  column,  if  not  an 
'error  will  be  written  and  the  user  will  be  notified. 

t 

'@param  -  CheckByOccField  is  a  boolean  variable  that  will  be  True  if  the  reenlistment 
rates  need  to  be  verified  by  OccField  or  false  if  the  rates  need  to  be 
verified  by  MOS 

t 

Sub  ReenlistmentRateCheckMaster(bCheckByOccField  As  Boolean) 

Dim  lLength  As  Long 
Dim  lWidth  As  Long 
'Data  Sheet  Ranges 
Dim  rZoneAMOS  As  Range 
Dim  rZoneAOccField  As  Range 
Dim  rZoneAEASPop  As  Range 
Dim  rZoneAB  Space  As  Range 
Dim  rZoneARRR  As  Range 
Dim  rZoneBMOS  As  Range 
Dim  rZoneBOccField  As  Range 
Dim  rZoneBEASPop  As  Range 
Dim  rZoneBBSpace  As  Range 
Dim  rZoneBRRR  As  Range 
Dim  rZoneCMOS  As  Range 
Dim  rZoneCOccField  As  Range 
Dim  rZoneCEASPop  As  Range 
Dim  rZoneCBSpace  As  Range 
Dim  rZoneCRRR  As  Range 
'Reenlistment  Rates  Ranges 
Dim  rRRRZoneAOccField  As  Range 
Dim  rRRRZoneBOccField  As  Range 
Dim  rRRRZoneCOccField  As  Range 
Dim  rRRRZoneA  As  Range 
Dim  rRRRZoneB  As  Range 
Dim  rRRRZoneC  As  Range 
'Check  Zone  A  by  OccField 

Set  rZoneAMOS  =  wsMOSZoneA.Range("ZoneAMOSStart",  wsMOSZoneA.Range 
("ZoneAMOSStart").End(xlDown)) 
lLength  =  rZoneAMOS. Rows. Count 

Set  rZoneAOccField  =  wsMOSZoneA.Range("ZoneAOccFieldStart",  _ 
wsMOSZoneA.Range("ZoneAOccFieldStart")  _ 
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.Offset(lLength,  0)) 

Set  rZoneAEASPop  =  wsMOSZoneA.Range("ZoneAEASPopStart",  _ 
wsMOSZoneA.Range("ZoneAEASPopStart")  _ 

.End(xlDown)) 

Set  rZoneABSpace  =  wsMOSZoneA.Range("ZoneABSAdjStart",  _ 
wsMOSZoneA.Range("ZoneABSAdj Start")  _ 

.End(xlDown)) 

Set  rRRRZoneAOccField  =  wsRRRZoneA.Range("RRRZoneAOccFieldStart",  _ 
wsRRRZoneA.Range("RRRZoneAOccFieldStart")  _ 

.End(xlDown)) 

Set  rRRRZoneA  =  wsRRRZoneA.Range(wsRRRZoneA.Range("RRRZoneAStart") 
.End(xlToRight),  wsRRRZoneA.Range("RRRZoneAStart")  _ 
.End(xlDown)) 

lWidth  =  rRRRZoneA. Columns. Count 

Set  rZoneARRR  =  wsMOSZoneA.Range("ZoneARRRStart",  wsMOSZoneA. Range 
("ZoneARRRStart").Offset(lLength  +  1,  lWidth+  1)) 

Call  CheckRRR(rZoneAOccField,  rZoneAMOS,  rZoneAEASPop,  rZoneABSpace,  _ 
rZoneARRR,  rRRRZoneAOccField,  rRRRZoneA,  lLength,  _ 
lWidth,  bCheckByOccField) 

If  INumErrors  >  0  Then 

MsgBox  ("Error:  There  are  OccField/MOS(s)  in  Zone  A  not  listed  in  the  "  _ 

&  "Reenlistment  Rates-Zone  A  worksheet."  &  vbCrLf  &  "  See  Data  -  "  _ 

&  "Zone  A  Workseet.") 

End  If 

'Check  Zone  B  by  OccField 

Set  rZoneBMOS  =  wsMOSZoneB.Range("ZoneBMOSStart",  wsMOSZoneB. Range 
("ZoneBMOSStart").End(xlDown)) 
lLength  =  rZoneBMOS. Rows. Count 

Set  rZoneBOccField  =  wsMOSZoneB. Range("ZoneBOccFieldStart",  _ 
wsMOSZoneB. Range("ZoneBOccFieldStart")  _ 

.Offset(lLength,  0)) 

Set  rZoneBEASPop  =  wsMOSZoneB. Range("ZoneBEASPopStart",  _ 
wsMOSZoneB.Range("ZoneBEASPopStart")  _ 

.End(xlDown)) 

Set  rZoneBBSpace  =  wsMOSZoneB. Range("ZoneBBSAdjStart",  _ 
wsMOSZoneB. Range("ZoneBBSAdj Start")  _ 

.End(xlDown)) 

Set  rRRRZoneBOccField  =  wsRRRZoneB.Range("RRRZoneBOccFieldStarf ',  _ 
wsRRRZoneB.Range("RRRZoneBOccFieldStart")  _ 

.End(xlDown)) 

Set  rRRRZoneB  =  wsRRRZoneB.Range(wsRRRZoneB.Range("RRRZoneBStarf ')  _ 
.End(xlToRight),  wsRRRZoneB.Range("RRRZoneBStart")  _ 
.End(xlDown)) 

lWidth  =  rRRRZoneB. Columns. Count 

Set  rZoneBRRR  =  wsMOSZoneB. Range("ZoneBRRRStart",  wsMOSZoneB. Range 
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("ZoneBRRRStart").Offset(lLength  +  1,  lWidth  +  1)) 

Call  CheckRRR(rZoneBOccField,  rZoneBMOS,  rZoneBEASPop,  rZoneBBSpace,  _ 
rZoneBRRR,  rRRRZoneBOccField,  rRRRZoneB,  lLength,  _ 
lWidth,  bCheckByOccField) 

If  INumErrors  >  0  Then 

MsgBox  ("Error:  There  are  OccField/MOS(s)  in  Zone  B  not  listed  in  the  "  _ 

&  "Reenlistment  Rates-Zone  B  worksheet."  &  vbCrLf  &  "  See  Data  -  "  _ 

&  "Zone  B  Workseet.") 

End  If 

'Check  Zone  C  by  OccField 

Set  rZoneCMOS  =  wsMOSZoneC.Range("ZoneCMOSStart",  wsMOSZoneC. Range 
("ZoneCMOSStarf').End(xlDown)) 
lLength  =  rZoneCMOS. Rows. Count 

Set  rZoneCOccField  =  wsMOSZoneC. Range("ZoneCOccFieldStart",  _ 
wsMOSZoneC. Range("ZoneCOccFieldStart")  _ 

.Offset(lLength,  0)) 

Set  rZoneCEASPop  =  wsMOSZoneC. Range("ZoneCEASPopStart",  _ 
wsMOSZoneC.Range("ZoneCEASPopStart")  _ 

.End(xlDown)) 

Set  rZoneCBSpace  =  wsMOSZoneC. Range("ZoneCBSAdjStart",  _ 
wsMOSZoneC. Range("ZoneCBSAdj Start")  _ 

.End(xlDown)) 

Set  rRRRZoneCOccField  =  wsRRRZoneC.Range("RRRZoneCOccFieldStarf _ 
wsRRRZoneC.Range("RRRZoneCOccFieldStart")  _ 

.End(xlDown)) 

Set  rRRRZoneC  =  wsRRRZoneC.Range(wsRRRZoneC.Range("RRRZoneCStarf ')  _ 
.End(xlToRight),  wsRRRZoneC.Range("RRRZoneCStart")  _ 
.End(xlDown)) 

lWidth  =  rRRRZoneC. Columns. Count 

Set  rZoneCRRR  =  wsMOSZoneC. Range("ZoneCRRRStart",  wsMOSZoneC. Range 
("ZoneCRRRStart").Offset(lLength  +  1,  lWidth  +  1)) 

Call  CheckRRR(rZoneCOccField,  rZoneCMOS,  rZoneCEASPop,  rZoneCBSpace,  _ 
rZoneCRRR,  rRRRZoneCOccField,  rRRRZoneC,  lLength,  _ 
lWidth,  bCheckByOccField) 

If  INumErrors  >  0  Then 

MsgBox  ("Error:  There  are  OccField/MOS(s)  in  Zone  C  not  listed  in  the  "  _ 

&  "Reenlistment  Rates-Zone  C  worksheet."  &  vbCrLf  &  "  See  Data  -  "  _ 

&  "Zone  C  Workseet.") 

End  If 
End  Sub 

'CheckRRRQ  is  a  subroutine  that  will  verify  that  the  MOS/OccField  is  listed  in  the 
'Reenlistment  Rate  Worksheet  for  the  correct  zone.  If  the  MOS/OccField  is  listed 
'the  reenlistment  rates  and  its  corresponding  multiplier  will  be  listed  in  the 
'reenlistment  rate  column, if  not  an  error  will  be  written. 
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'@param  -  rOccFieldZoneList  is  the  range  where  the  correct  OccField  will  be  listed 
on  the  Data  Sheet 

-  rZoneMOSList  is  the  range  of  MOS's  that  need  to  be  verified 

-  rEASPop  is  the  range  of  EAS  Populations  corresponding  to  the  MOS  that  is 

'  to  be  verified 

-  rBSpace  is  the  range  of  Boat  Spaces  available  for  the  MOS  that  is  to  be 

'  verified 

"  -  rZoneRRR  is  the  range  that  will  be  filled  with  reenlistment  rates  from 

'  the  Reenlistment  Rate  Worksheet 

-  rRRRZoneList  is  the  range  of  MOS  or  OccFields  that  the  ZoneMOSList  will  be 
verified  against 

'  -  rRRRZone  is  the  range  of  reenlistment  rates  corresponding  to  the  RRRZoneList 

-  INumToCheck  is  the  number  of  MOS/OccFields  to  check 

-  INumMult  is  the  number  of  SRB  multipliers  that  are  being  used 

-  bCheckBy  OccField  is  a  boolean  variable  that  will  be  True  if  the  reenlistment 
rates  need  to  be  verified  by  OccField  or  false  if  the  rates  need  to  be 

'  verified  by  MOS 

t 

Sub  CheckRRR(rOccFieldZoneList  As  Range,  rZoneMOSList  As  Range,  _ 
rEASPop  As  Range,  rBSpace  As  Range,  rZoneRRR  As  Range,  _ 
rRRRZoneList  As  Range,  rRRRZone  As  Range,  INumToCheck  As  Long,  _ 
INumMult  As  Long,  bCheckByOccField  As  Boolean) 

Dim  INumRRR  As  Long 
Dim  lPoint  As  Long 
Dim  sOccField  As  String 
Dim  bContain  As  Boolean 
rOccF  ieldZoneList.  Clear 
rZoneRRR.  Clear 
INumErrors  =  0 
INumNodes  =  0 

INumRRR  =  rRRRZoneList.Rows. Count 

IMaxIndex  =  3  *  INumRRR 

ReDim  llndex(l  To  IMaxIndex)  As  Long 

ReDim  sNodeNames(0  To  2  *  INumRRR)  As  String 

Dim  i  As  Integer 

Dim  j  As  Integer 

Dim  lRRRPointer()  As  Long 

ReDim  lRRRPointer(  1  To  IMaxIndex) 

For  i  =  1  To  INumRRR 
AddNode  rRRRZoneList(i) 

lRRRPointer(LookupNode(rRRRZoneList(i)))  =  i  +  1 
Next  i 

Formating  here  so  the  0's  on  the  MOSs  will  show  up 
Call  FormatOccFieldText(rOccFieldZoneList) 
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'Prints  the  multiplier  Numbers 
For  i  =  1  To  INumMult 

rZoneRRR(  1 ,  i)  =  rRRRZone(  1 ,  i) 

Next  i 

For  i  =  1  To  INumToCheck 

sOccField  =  Mid(rZoneMOSList(i),  1,  2) 

If  bCheckByOccField  Then 

bContain  =  DoesContain(sOccField) 

If  bContain  Then 

lPoint  =  lRRRPointer(LookupNode(sOccField)) 

End  If 
Else 

bContain  =  DoesContain(rZoneMOSList(i)) 

If  bContain  Then 

lPoint  =  lRRRPointer(LookupNode(rZoneMOSList(i))) 

End  If 
End  If 

rOccFieldZoneList(i)  =  sOccField 
If  bContain  Then 

For  j  =  1  To  INumMult 

rZoneRRR(i  +  1,  j)  =  Application. WorksheetFunction.Round  _ 
(rRRRZone(lPoint,  j),  4) 

Next  j 
Else 

If  rEASPop(i)  =  0  Then 

rZoneRRR(i  +1,1)  =  "Error,  MOS  not  LISTED  and  NOT  needed  due  " 
&  "to  Zero  EAS  Population" 

Elself  rBSpace(i)  =  0  Then 

rZoneRRR(i  +  1,  1)  =  "Error,  MOS  not  LISTED  and  NOT  needed  due  " 
&  "to  Zero  Boat  Spaces" 

Else 

rZoneRRR(i  +1,1)  =  "Error,  MOS  not  LISTED  and  NEEDED" 
INumErrors  =  INumErrors  +  1 
End  If 
End  If 
Next  i 

Call  FormatTwoDecimal(rZoneRRR) 

End  Sub 

'File:  HashCode  Module  in  VBA 
'Date  Created:  March  2007 
'Last  Updated:  May  2007 

t 

'Checks  is  a  VBA  Module  that  contains  a  list  of  subroutines  that  will  perform 
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'checks  on  the  data  set  to  verify  that  it  is  complete  in  order  to  run  the 
'optimization  software. 

t 

Author  Matthew  Catyle 
'  Updated  by  Kent  Robbins,  Jr. 

t 

Option  Explicit 

'HashName  is  a  function  that  will  take  a  string  variable  and  convert  it  to  a 
'numerical  value  for  use  in  a  hash  code. 

f 

'@param  -  sName  is  a  String  Variable  that  will  be  converted  to  a  numerical  value 
'@return  -  HashName  is  a  Long  variable  that  is  converted  from  sName 

t 

Function  HashName(sName  As  String)  As  Long 

'Converts  an  input  string  sName  into  a  long  integer  between  1  and  IMaxIndex 
Dim  lValue  As  Long 
Dim  i  As  Integer 
lValue  =  0 

For  i  =  1  To  Len( sName) 

lValue  =  (10  *  lValue  +  Asc(Mid( sName,  i,  1)))  Mod  IMaxIndex 
Next  i 

HashName  =  lValue  +  1 
End  Function 

'AddNode  is  a  subroutine  that  will  add  a  variable  to  the  hash  code. 

t 

'@param  -  sName  is  a  string  variable  that  will  be  added  to  the  hash  code. 

t 

Sub  AddNode(sName  As  String) 

Dim  hv  As  Long 

hv  =  HashName(sName) 

Do  While  llndex(hv)  <>  0  And  sNodeNames(lIndex(hv))  <>  sName 
hv  =  hv  +  1 

If  hv  >  IMaxIndex  Then 
hv  =  1 
End  If 
Loop 

If  llndex(hv)  =  0  Then 

INumNodes  =  INumNodes  +  1 
llndex(hv)  =  INumNodes 
sNodeNames(lNumNodes)  =  sName 
End  If 
End  Sub 

'LoopupNode  is  a  function  that  will  return  the  node  number  of  the  parameter. 
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'@param  -  sName  is  a  string  variable  used  to  look  up  its  associated  node  number. 
'@return  -  LoopUpNode 

t 

Function  LookupNode(sName  As  String)  As  Long 
Dim  hv  As  Long 
hv  =  HashName(  sName) 

Do  While  llndex(hv)  <>  0  And  sNodeNames(lIndex(hv))  <>  sName 
hv  =  hv  +  1 

If  hv  >  IMaxIndex  Then 
hv  =  1 
End  If 
Loop 

LookupNode  =  llndex(hv) 

End  Function 

'DoesContain  is  a  function  that  will  compute  the  correct  number  of  reenlistments  for 
'Zone  A  that  are  needed  when  BNA  school  seats  are  used. 

t 

'@param  -  sName  is  a  string  variable  that  needs  to  be  check  to  see  if  it  is  in  the 
hash  code. 

'@return  -  DoesContain  is  a  boolean  variable.  True  if  sName  is  in  the  hash  code, 
otherwise  false. 

t 

Function  DoesContain( sName  As  String)  As  Boolean 
Dim  hv  As  Long 
hv  =  HashName(sName) 

Do  While  llndex(hv)  <>  0  And  sNodeNames(lIndex(hv))  <>  sName 
hv  =  hv  +  1 

If  hv  >  IMaxIndex  Then 
hv  =  1 
End  If 
Loop 

If  llndex(hv)  =  0  Then 
DoesContain  =  False 
Else 

DoesContain  =  True 
End  If 

End  Function 

'File:  Reduction  Module  in  VBA 
'Date  Created:  1 7  June  2007 
'Last  Updated:  25  July  2007 

t 

'Reduction  is  a  VBA  Modula  that  will  combine  all  MOS  from  Zone  A,  B,  and  C  onto  a 
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'single  worksheet  containing  the  necessary  data  for  the  GAMS  and  VBA  solver  to  find 
'the  Optimum  SRB  multiplier.  Reduction  will  only  transfer  MOSs  that  will  not  achieve 
'the  reenlistment  goal  without  the  use  of  a  multiplier  greater  than  0.  Reduction 
'will  also  check  the  cost  of  each  multiplier  and  will  assign  it  a  0.0  reenlistment  rate 
'if  that  multiplier  violates  the  maximum  bonus  contstraint. 

t 

'@  Author  Kent  Robbins,  Jr. 

t 

Option  Explicit 

'ReduceSizeMaster  is  a  subroutine  that  will  gather  the  data  neccessary  and  call  the 
function  Reduce  to  execute  the  reduction  and  transfer  of  the  MOSsin  Zones  A,  B,  and 
'C  that  will  not  achieve  the  reenlistment  goal  with  a  multiplier  of  0.  Results  are 
'transfered  to  the  ReducedFull  Problem  Worksheet. 

t 

'@param  -  none 

t 

Sub  ReduceSizeMaster() 

'Zone  Data  Ranges 
Dim  rZoneAOccField  As  Range 
Dim  rZoneAMOSField  As  Range 
Dim  rZoneAZoneField  As  Range 
Dim  rZoneAMOSPopField  As  Range 
Dim  rZoneAEASPopField  As  Range 
Dim  rZoneAPreMultField  As  Range 
Dim  rZoneAWeightField  As  Range 
Dim  rZoneAReupNeedField  As  Range 
Dim  rZoneATrainCostField  As  Range 
Dim  rZoneARRR  As  Range 
Dim  dZoneAAvgPay  As  Double 
Dim  rZoneBOccField  As  Range 
Dim  rZoneBMOSField  As  Range 
Dim  rZoneBZoneField  As  Range 
Dim  rZoneBMOSPopField  As  Range 
Dim  rZoneBEASPopField  As  Range 
Dim  rZoneBPreMultField  As  Range 
Dim  rZoneBWeightField  As  Range 
Dim  rZoneBReupNeedField  As  Range 
Dim  rZoneBTrainCostField  As  Range 
Dim  rZoneBRRR  As  Range 
Dim  dZoneBAvgPay  As  Double 
Dim  rZoneCOccField  As  Range 
Dim  rZoneCMOSField  As  Range 
Dim  rZoneCZoneField  As  Range 
Dim  rZoneCMOSPopField  As  Range 
Dim  rZoneCEASPopField  As  Range 
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Dim  rZoneCPreMultField  As  Range 
Dim  rZoneCWeightField  As  Range 
Dim  rZoneCReupNeedField  As  Range 
Dim  rZoneCTrainCostField  As  Range 
Dim  rZoneCRRR  As  Range 
Dim  dZoneCAvgPay  As  Double 
Dim  dMaxIndivBonus  As  Double 
'OutPutZone  Variables 
Dim  iMaxOutputLength  As  Integer 
Dim  iMaxOutputWidth  As  Integer 
Dim  rOutput  As  Range 
Dim  rAvgPays  As  Range 
Dim  i  As  Integer 
Dim  iLinePosition  As  Integer 

dMaxIndivBonus  =  wsReducedProblem.Range("MaxIndBonus") 

Set  rZoneAOccField  =  wsMOSZoneA.Range("ZoneAOccFieldStart",  _ 
wsMOSZoneA.Range("ZoneAOccFieldStart")  _ 

.End(xlDown)) 

Set  rZoneAMOSField  =  wsMOSZoneA.Range("ZoneAMOSStart",  _ 
wsMOSZoneA.Range("ZoneAMOSStart")  _ 

.End(xlDown)) 

Set  rZoneAZoneField  =  wsMOSZoneA.Range("ZoneASRBZoneStart",  _ 
wsMOSZoneA.Range("ZoneASRBZoneStart")  _ 

.End(xlDown)) 

Set  rZoneAMOSPopField  =  wsMOSZoneA.Range("ZoneAMOSPopStart",  _ 
wsMOSZoneA.Range("ZoneAMOSPopStart")  _ 

.End(xlDown)) 

Set  rZoneAEASPopField  =  wsMOSZoneA.Range("ZoneAEASPopStart",  _ 
wsMOSZoneA.Range("ZoneAEASPopStart")  _ 

.End(xlDown)) 

Set  rZoneAPreMultField  =  wsMOSZoneA.Range("ZoneAPresetStart",  _ 
wsMOSZoneA.Range("ZoneAPresetStart")  _ 

.End(xlDown)) 

Set  rZoneAWeightField  =  wsMOSZoneA.Range("ZoneAWgtFctStart",  _ 
wsMOSZoneA.Range("ZoneAWgtFctStart")  _ 

.End(xlDown)) 

Set  rZoneAReupNeedField  =  wsMOSZoneA.Range("ZoneABSAdj Start",  _ 
wsMOSZoneA.Range("ZoneABS Adj Start")  _ 

.End(xlDown)) 

Set  rZoneATrainCostField  =  wsMOSZoneA.Range("ZoneATmCstStart",  _ 
wsMOSZoneA.Range("ZoneATmCstStart")  _ 

.End(xlDown)) 

Set  rZoneARRR  =  wsMOSZoneA.Range(wsMOSZoneA.Range("ZoneARRRStart") 
.End(xlToRight),  wsMOSZoneA.Range("ZoneARRRStart")  _ 
.End(xlDown)) 
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dZoneAAvgPay  =  wsMOSZoneA.Range("ZoneAAvgPay") 

Set  rZoneBOccField  =  wsMOSZoneB.Range("ZoneBOccFieldStart",  _ 
wsMOSZoneB.Range("ZoneBOccFieldStart")  _ 

.End(xlDown)) 

Set  rZoneBMOSField  =  wsMOSZoneB.Range("ZoneBMOSStart",  _ 
wsMOSZoneB.Range("ZoneBMOSStart")  _ 

.End(xlDown)) 

Set  rZoneBZoneField  =  wsMOSZoneB.Range("ZoneBSRBZoneStart",  _ 
wsMOSZoneB.Range("ZoneBSRBZoneStart")  _ 

.End(xlDown)) 

Set  rZoneBMOSPopField  =  wsMOSZoneB.Range("ZoneBMOSPopStart",  _ 
wsMOSZoneB  .Range("ZoneBMOSPopStart")  _ 

.End(xlDown)) 

Set  rZoneBEASPopField  =  wsMOSZoneB. Range("ZoneBEASPopStart",  _ 
wsMOSZoneB  .Range("ZoneBEASPopStart")  _ 

.End(xlDown)) 

Set  rZoneBPreMultField  =  wsMOSZoneB.Range("ZoneBPresetStart",  _ 
wsMOSZoneB. Range("ZoneBPresetStart")  _ 

.End(xlDown)) 

Set  rZoneBWeightField  =  wsMOSZoneB. Range("ZoneBWgtFctStart",  _ 
wsMOSZoneB. Range("ZoneBWgtFctStart")  _ 

.End(xlDown)) 

Set  rZoneBReupNeedField  =  wsMOSZoneB.Range("ZoneBBSAdjStart",  _ 
wsMOSZoneB. Range("ZoneBBSAdjStart")  _ 

.End(xlDown)) 

Set  rZoneBTrainCostField  =  wsMOSZoneB.Range("ZoneBTrnCstStart",  _ 
wsMOSZoneB. Range("ZoneBTrnCstStart")  _ 

.End(xlDown)) 

Set  rZoneBRRR  =  wsMOSZoneB. Range(wsMOSZoneB.Range("ZoneBRRRStart") 
.End(xlToRight),  wsMOSZoneB. Range("ZoneBRRRStart")  _ 
.End(xlDown)) 

dZoneBAvgPay  =  wsMOSZoneB. Range("ZoneBAvgPay") 

Set  rZoneCOccField  =  wsMOSZoneC.Range("ZoneCOccFieldStart",  _ 
wsMOSZoneC.Range("ZoneCOccFieldStart")  _ 

.End(xlDown)) 

Set  rZoneCMOSField  =  wsMOSZoneC.Range("ZoneCMOSStart",  _ 
wsMOSZoneC.Range("ZoneCMOSStart")  _ 

.End(xlDown)) 

Set  rZoneCZoneField  =  wsMOSZoneC.Range("ZoneCSRBZoneStart",  _ 
wsMOSZoneC.Range("ZoneCSRBZoneStart")  _ 

.End(xlDown)) 

Set  rZoneCMOSPopField  =  wsMOSZoneC.Range("ZoneCMOSPopStart",  _ 
wsMOSZoneC.Range("ZoneCMOSPopStart")  _ 

.End(xlDown)) 

Set  rZoneCEASPopField  =  wsMOSZoneC.Range("ZoneCEASPopStart",  _ 
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wsMOSZoneC.Range("ZoneCEASPopStart")  _ 

.End(xlDown)) 

Set  rZoneCPreMultField  =  wsMOSZoneC.Range("ZoneCPresetStart",  _ 
wsMOSZoneC.Range("ZoneCPresetStart")  _ 

.End(xlDown)) 

Set  rZoneCWeightField  =  wsMOSZoneC.Range("ZoneCWgtFctStart",  _ 
wsMOSZoneC.Range("ZoneCWgtFctStart")  _ 

.End(xlDown)) 

Set  rZoneCReupNeedField  =  wsMOSZoneC.Range("ZoneCBSAdjStart",  _ 
wsMOSZoneC.Range("ZoneCBSAdjStart")  _ 

.End(xlDown)) 

Set  rZoneCTrainCostField  =  wsMOSZoneC.Range("ZoneCTrnCstStart",  _ 
wsMOSZoneC.Range("ZoneCTrnCstStart")  _ 

.End(xlDown)) 

Set  rZoneCRRR  =  wsMOSZoneC.Range(wsMOSZoneC.Range("ZoneCRRRStart") 
.End(xlToRight),  wsMOSZoneC.Range("ZoneCRRRStart")  _ 
.End(xlDown)) 

dZoneCAvgPay  =  wsMOSZoneC.Range("ZoneCAvgPay") 

iMaxOutputLength  =  rZoneAOccField.Rows.  Count  +  rZoneBOccField.Rows.  Count 
+  rZoneCOccField.Rows. Count 
iMaxOutputWidth  =  9  +  rZoneARRR.Coluinns. Count 
iLinePosition  =  2 

Set  rOutput  =  wsReducedProblem.Range("RFPOutputStart",  _ 
wsReducedProblem.Range("RFPOutputStart")  _ 

.Offset(iMaxOutputLength  +  1,  iMaxOutputWidth)) 

Set  rAvgPays  =  wsReducedProblem.Range("AvgPayA",  wsReducedProblem.Range 
("AvgPayA").Offset(3,  0)) 
rOutput.  Clear 

Call  FormatOccFieldText(rOutput) 
rOutput.HorizontalAlignment  =  xlCenter 
rAvgPays. Clear 

rAvgPays.HorizontalAlignment  =  xlCenter 

rOutput(l,  1)  =  "Occ  Field" 

rOutput(  1 ,  2)  =  "MOS" 

rOutput(l,  3)  =  "Zone" 

rOutput(  1 ,  4)  =  "MOSPop" 

rOutput(  1 ,  5)  =  "EASPop" 

rOutput(  1,6)  =  "PresetMult" 

rOutput(l,  7)  =  "WgtFactor" 

rOutput(  1 ,  8)  =  "B/SPACE" 

rOutput(l,  9)  =  "TrnCosts" 

For  i  =  1  To  rZoneARRR.Columns. Count 
rOutput(  1 ,  9  +  i)  =  rZoneARRR(  1 ,  i) 

Next  i 

rAvgPays(l)  =  dZoneAAvgPay 
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rAvgPays(2)  =  dZoneBAvgPay 
rAvgPays(3)  =  dZoneCAvgPay 

iLinePosition  =  Reduce(rZoneAOccField,  rZoneAMOSField,  rZoneAZoneField,  _ 
rZoneAMOSPopField,  rZoneAEASPopField,  rZoneAPreMultField,  _ 
rZoneAWeightField,  rZoneAReupNeedField,  rZoneATrainCostField,  _ 
rZoneARRR,  rOutput,  dZoneAAvgPay,  dMaxIndivBonus,  iLinePosition) 
iLinePosition  =  Reduce(rZoneBOccField,  rZoneBMOSField,  rZoneBZoneField,  _ 
rZoneBMOSPopField,  rZoneBEASPopField,  rZoneBPreMultField,  _ 
rZoneBWeightField,  rZoneBReupNeedField,  rZoneBTrainCostField,  _ 
rZoneBRRR,  rOutput,  dZoneBAvgPay,  dMaxIndivBonus,  iLinePosition) 
iLinePosition  =  Reduce(rZoneCOccField,  rZoneCMOSField,  rZoneCZoneField,  _ 
rZoneCMOSPopField,  rZoneCEASPopField,  rZoneCPreMultField,  _ 
rZoneCWeightField,  rZoneCReupNeedField,  rZoneCTrainCostField,  _ 
rZoneCRRR,  rOutput,  dZoneCAvgPay,  dMaxIndivBonus,  iLinePosition) 

End  Sub 

'Reduce  is  a  function  that  will  execute  the  reduction  and  transfer  of  those  MOSs  in 
'a  Zone  that  will  need  a  multiplier  greater  than  0  to  possible  achieve  the  reenlistment 
'goal.  Reduce  will  also  check  the  cost  of  using  each  multiplier  and  if  the  cost  is 
'greater  than  the  maximum  allowable  bonus  the  reenlistment  rate  for  that  multiplier  will 
'be  set  to  0  so  it  will  not  be  used. 

t 

'@param  -  rOccField  is  the  range  of  OccFields  of  the  Zones  MOSs  to  be  reduced. 

'  -  rMos  is  the  range  of  MOSs  in  the  Zone  that  is  to  reduced. 

-  rZone  is  a  range  variable  and  is  the  Zone  that  to  be  reduced  for  each  MOS. 

-  rMosPop  is  a  range  variable  of  MOS  populations  for  each  MOS. 

-  rEASPop  is  a  range  variable  of  EAS  populations  for  each  MOS. 

-  rPreset  is  a  range  variable  of preset  multiplier  value  for  each  MOS. 

-  rWeight  is  a  range  variable  of  exogenous  weighting  factor  for  each  MOS. 

-  rBSpace  is  a  range  variable  of  needed  reenlistment  numbers  for  each  MOS. 

-  rTrnCst  is  a  range  variable  of  training  costs  for  each  MOS. 

'  -  rRRR  is  a  range  variable  of  reenlistment  response  rates  for  each  multiplier. 

-  rOutput  is  a  range  where  the  MOS  and  data  will  be  transfered  to. 

-  dAvgPay  is  a  double  variable  of  the  average  pay  of  the  Zone  being  reduced. 

’  -  dMaxBonus  is  a  double  variable  of  the  maximum  allowable  bonus. 

'  -  iLine  is  an  integer  variable  and  a  place  holder  for  the  location  in  the 

output  for  the  next  Zone  to  start  being  transfered  to. 

t 

'@return  -  Reduce  is  an  integer  variable  and  is  the  location  marker  for  the  output 
'  for  the  next  Zone  that  is  to  be  reduced. 

t 

Function  Reduce(rOccField  As  Range,  rMos  As  Range,  rZone  As  Range,  _ 
rMosPop  As  Range,  rEASPop  As  Range,  rPreset  As  Range,  _ 
rWeight  As  Range,  rBSpace  As  Range,  rTrnCst  As  Range,  _ 
rRRR  As  Range,  rOutput  As  Range,  dAvgPay  As  Double,  _ 
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dMaxBonus  As  Double,  iLine  As  Integer)  As  Integer 
Dim  dGap  As  Double 
Dim  i  As  Integer 
Dim  j  As  Integer 
Dim  bContinue  As  Boolean 
dGap  =  rRRR(l,  2)  -  rRRR(  1 ,  1) 

For  i  =  1  To  rOccField.Rows. Count 
If  rEASPop(i)  =  0  Then 
bContinue  =  False 
Elself  rBSpace(i)  =  0  Then 
bContinue  =  False 
Else 

bContinue  =  True 
End  If 

If  bContinue  Then 
If  rPreset(i)  >  -1  Then 

rOutput(  iLine,  1)  =  rOccField(i) 
rOutput(iLine,  2)  =  rMos(i) 
rOutput(iLine,  3)  =  rZone(i) 
rOutput(iLine,  4)  =  rMosPop(i) 
rOutput(iLine,  5)  =  rEASPop(i) 
rOutput(iLine,  6)  =  rPreset(i) 
rOutput( iLine,  7)  =  rWeight(i) 
rOutput(iLine,  8)  =  rBSpace(i) 
rOutput(iLine,  9)  =  rTrnCst(i) 

For  j  =  1  To  rRRR.Columns. Count 
rOutput(iLine,  9  +  j)  =  rRRR(i  +  1,  j) 

Next  j 

iLine  =  iLine  +  1 

Elself  (rEASPop(i)  *  rRRR(i  +  1,  1)  <  rBSpace(i))  Then 
rOutput(  iLine,  1)  =  rOccField(i) 
rOutput(iLine,  2)  =  rMos(i) 
rOutput(iLine,  3)  =  rZone(i) 
rOutput(iLine,  4)  =  rMosPop(i) 
rOutput(iLine,  5)  =  rEASPop(i) 
rOutput(iLine,  6)  =  rPreset(i) 
rOutput( iLine,  7)  =  rWeight(i) 
rOutput(iLine,  8)  =  rBSpace(i) 
rOutput(iLine,  9)  =  rTrnCst(i) 

For  j  =  1  To  rRRR.Columns. Count 

If  (dGap  *  (j  -  1)  *  dAvgPay  *  4  >  dMaxBonus)  Then 
rOutput(iLine,  9  +  j)  =  0 
Else 

rOutput(iLine,  9  +  j)  =  rRRR(i  +  1,  j) 

End  If 
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Next  j 

iLine  =  iLine  +  1 
End  If 
End  If 
Next  i 

Reduce  =  iLine 
End  Function 

'File:  So/veVBA  Module  in  VBA 
'Date  Created:  1 7  June  2007 
'Last  Updated:  31  August  2007 

t 

'Solve  VBA  Module  is  a  VBA  module  that  will  solve  the  Lagrangian  relaxation  problem  of 
'finging  the  optimum  SRB  multipliers  to  assign  to  each  MOSZ  based  on  the  weighted 
'penalty  function  occured  with  each  multiplier  and  subject  to  certain  buget 
'limitations. 

f 

'@  Author  Kent  Robbins,  Jr. 

t 

Option  Explicit 
Public  dlnfinity  As  Double 
Public  dBudget  As  Double 
Public  dTest  As  Double 
Public  iNumMOS  As  Integer 
Public  iNumMult  As  Integer 
Public  dGap  As  Double 
Public  dHLambda  As  Double 
Public  dZL  As  Double 
Public  dZU  As  Double 
Public  dZH  As  Double 
Public  dSumCost  As  Double 
Public  dFinalCost  As  Double 
Public  dTotalDeviation  As  Double 

'SolveEXCEL  is  the  master  subroutine  that  run  the  VBA  Model.  This  subroutine  will 
'gather  the  pertinent  data  that  is  need  in  later  subroutines  to  find  the  optimum 
'SRB  multipliers. 

t 

'@param  -  none 

t 

Sub  SolveEXCEL() 

Dim  dStartTime  As  Double 
dStartTime  =  Timer 
Dim  dFinishTime  As  Double 
'Worksheet  variables 
Dim  rMOSData  As  Range 
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Dim  rBudgetlnfo  As  Range 
Dim  dCost()  As  Double 
Dim  dDeviation()  As  Double 
Dim  i  As  Integer 
Dim  dSoln()  As  Double 
Dim  dPreset()  As  Double 
Dim  dOptimality  As  Double 
dTotalDeviation  =  0# 

'Read  in  data  from  worksheets 

Set  rMOSData  =  wsReducedProblem.Range(wsReducedProblem.Range  _ 
("RFPOutputStart").End(xlToRight),  wsReducedProblem.Range  _ 
("RFPOutputStart").End(xlDown)) 

Set  rBudgetlnfo  =  wsReducedProblem.Range("OverageWeight",  _ 

wsReducedProblem.Range("OverageWeight").End(xlDown)) 
dBudget  =  rBudgetInfo(2) 
dlnfinity  =  1  *  (10  A  22) 
iNumMOS  =  rMOSData.Rows. Count 
iNumMult  =  rMOSData.Columns. Count  -  9 
dGap  =  rMOSData(l,  1 1)  -  rMOSData(l,  10) 

ReDim  dCost(l  To  (iNumMOS  -  1),  1  To  iNumMult)  As  Double 
ReDim  dDeviation(l  To  (iNumMOS  -  1),  1  To  iNumMult)  As  Double 
ReDim  dPreset(l  To  (iNumMOS  -  1),  1  To  2)  As  Double 
ReDim  dSoln(l  To  (iNumMOS  -  1))  As  Double 
For  i  =  2  To  iNumMOS 

If  rMOSData(i,  6)  >=  0  Then 

dPreset(i  -  1,  1)  =  rMOSData(i,  6) 
dPreset(i  -  1,2)  =  rMOSData(i,  6)  /  dGap  +  1 
Else 

dPreset(i  -  1,  1)  =  -1 
dPreset(i  -  1 ,  2)  =  - 1 
End  If 
Next  i 

Call  FillCostandDeviation(rMOSData,  rBudgetlnfo,  dCost(),  dDeviation()) 

Call  Bound(dSoln(),  dPreset(),  dCost(),  dDeviation()) 

Call  Heuris(dDeviation(),  dCost(),  dSoln(),  dPreset()) 

Call  WriteResults(rMOSData,  rBudgetlnfo,  dSoln()) 

dFinishTime  =  Timer 

dOptimality  =  (dZH  -  dZL)  /  dZL  *  100 

MsgBox  ("VBA  Solve  has  finished."  &  vbCrLf  &  "Solution  is  within  "  &  Round 
(dOptimality,  6)  &  "%  of  optimality"  &  vbCrLf  &  "Run  time  was  "  &  _ 
Round((dFinishTime  -  dStartTime),  3)  &  "  seconds"  &  vbCrLf  &  "Object  " 
&  "Function  is  "  &  Round(dTotalDeviation,  4)  &  vbCrFf  &  "Total  Cost  "  _ 
&  "is  $"  &  Round(dFinalCost  *  1000,  2)) 

End  Sub 

'FiUCostandDeviation  is  a  subroutine  that  will  fill  the  Cost  and  Deviation  arrays 
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for  each  MOSZ  and  multiplier  combination. 


'@param  -  rMOSData  is  a  range  variable  that  contains  data  for  all  the  MOSZs. 

-  rBudgetlnfo  is  a  range  variable  that  contains  the  budget  and  pay 
information  needed  for  the  SRB  Optimization  problem. 

-  dCost  is  an  array  of  doubles  that  will  contain  the  cost  for  setting  each 
MOSZ  with  all  available  multipliers. 

-  dDeviation  is  an  array  of  doubles  that  will  contain  the  deviation  for 
setting  each  MOSZ  with  all  available  multipliers. 

t 

Sub  FillCostandDeviation(rMOSData  As  Range,  rBudgetlnfo  As  Range,  _ 
dCost()  As  Double,  dDeviation()  As  Double) 

Dim  i  As  Integer 
Dim  dScaleFactor  As  Double 
Dim  dWeight  As  Double 
Dim  dEligible  As  Double 
Dim  dNeededNum  As  Double 
Dim  dTotalNum  As  Double 
Dim  dTrainCost  As  Double 
Dim  dAvgPay  As  Double 
Dim  iContractLength  As  Integer 
Dim  dMaxTrainCost  As  Double 
dScaleFactor  =  1000# 
dMaxTrainCost  =  0# 
dBudget  =  dBudget  /  dScaleFactor 
dTest  =  dlnfinity  /  1 . 1 
dHLambda  =  -dlnfinity 
Set  max  Training  Cost 
For  i  =  2  To  iNumMOS 

dMaxTrainCost  =  Application.WorksheetFunction.Max(dMaxTrainCost,  _ 
rMOSData(i,  9)) 

Next  i 

For  i  =  2  To  iNumMOS 

Dim  dTotalWeight  As  Double 

Dim  j  As  Integer 

Dim  dDev  As  Double 

Dim  dHL  As  Double 

dTotalNum  =  rMOSData(i,  4) 

dEligible  =  rMOSData(i,  5) 

dWeight  =  rMOSData(i,  7) 

dNeededNum  =  rMOSData(i,  8) 

dTrainCost  =  rMOSData(i,  9) 

dAvgPay  =  rBudgetInfo(rMOSData(i,  3)  +  4) 

iContractLength  =  rBudgetInfo(3) 

If  (dTrainCost  <  1#)  Then 
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dTrainCost  =  50# 

End  If 

If  (dTotalNum  >  0#)  Then 

dTotalWeight  =  dWeight  *  (dTrainCost  /  dMaxTrainCost)  /  dTotalNum 
Else 

dTotalWeight  =  dWeight  *  (dTrainCost  /  dMaxTrainCost) 

End  If 

If  (dEligible  <  1#)  Then 
dEligible  =  1# 

End  If 

For  j  =  1  To  iNumMult 

dDev  =  dNeededNum  -  dEligible  *  rMOSData(i,  9  +  j) 

If  (dDev  <  0#)  Then 

dDeviation(i  -  1,  j)  =  dTotalWeight  *  rBudgetlnfo(l)  *  (dDev  A  2) 

Else 

dDeviation(i  -  1,  j)  =  dTotalWeight  *  (dDev  A  2) 

End  If 

If  (rMOSData(i,  9  +  j)  =  0)  Then 
dCost(i  -  1,  j)  =  dlnfinity 
Else 

dCost(i  -  1,  j)  =  ((dEligible  *  rMOSData(i,  9  +  j))  *  dAvgPay  *  _ 
iContractLength  *  rMOSData(l,  9  +  j))  /  _ 
dScaleFactor 

End  If 

If  (j  >  1)  Then 

dHL  =  (dDeviation(i  -  1,  1)  -  dDeviation(i  -  1,  j))  /  dCost(i  -  1,  j) 

If  (dHL  >  dHLambda)  Then 
dHLambda  =  dHL 
End  If 
End  If 
Next  j 
Next  i 
End  Sub 

’Bound  is  a  subroutine  that  will  determine  the  upper  and  lower  bounds  of  the 
'Lagrangian  relaxation  and  call  two  other  subroutines  to  find  feasible  solution 
'sets  with  given  Lagrangian  multipliers  and  find  the  best  feasible  solution  set 
'that  has  the  lowest  object  function. 

t 

'@param  -  dSoln  is  an  array  of  doubles  that  contains  the  multipliers  solution  set. 
'  -  dPreset  is  an  array  of  doubles  the  list  the  preset  multipliers  if  chosen. 

-  dCost  is  an  array  of  doubles  that  will  contain  the  cost  for  setting  each 
MOSZ  with  all  available  multipliers. 

-  dDeviation  is  an  array  of  doubles  that  will  contain  the  deviation  for 
settingeach  MOSZ  with  all  available  multipliers. 
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Sub  Bound(dSoln()  As  Double,  dPreset()  As  Double,  dCost()  As  Double, 
dDeviation()  As  Double) 

Dim  dEps  As  Double 
Dim  dEndL  As  Double 
Dim  dEndR  As  Double 
Dim  dZUBest  As  Double 
Dim  dAmbda  As  Double 
Dim  dFLambda  As  Double 
Dim  i  As  Integer 
dEps  =  dHLambda  /  100000 
dEndL  =  0# 

dEndR  =  0.01  *  dHLambda 
dZUBest  =  dlnfinity 
i  =  0 
Do 

dAmbda  =  (dEndL  +  dEndR)  /  2# 

Call  MinFcn(dAmbda,  dPreset(),  dCost(),  dDeviation()) 

If  (dSumCost  <=  dBudget)  Then 
dEndR  =  dAmbda 
If  (dZU  <=  dZUBest)  Then 
dZUBest  =  dZU 
dFLambda  =  dAmbda 
End  If 
Exit  Do 
End  If 

dEndL  =  dAmbda 
dEndR  =  10.01  *  dEndR 
i  =  i  +  1 

Loop  Until  (i  =  2) 

Do 

dAmbda  =  (dEndL  +  dEndR)  /  2# 

Call  MinFcn( dAmbda,  dPreset(),  dCost(),  dDeviation()) 

If  (dSumCost  <=  dBudget)  Then 
dEndR  =  dAmbda 
If  (dZU  <=  dZUBest)  Then 
dZUBest  =  dZU 
dFLambda  =  dAmbda 
End  If 
Else 

dEndL  =  dAmbda 
End  If 

Loop  Until  ((dEndR  -  dEndL)  <=  dEps) 

Call  MinFeasible(dFLambda,  dSoln(),  dPreset(),  dDeviation(),  dCost()) 
End  Sub 
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'MinFcn  is  a  subroutine  that  will  determine  the  minimum  objective  function  for 
'the  current  Lagrangian  multiplier. 


'@param  -  dAmbda  is  the  Lagrangian  multiplier  current  used. 

-  dPreset  is  an  array  of  doubles  the  list  the  preset  multipliers  if  chosen. 

-  dCost  is  an  array  of  doubles  that  will  contain  the  cost  for  setting  each 
MOSZ  with  all  available  multipliers. 

-  dDeviation  is  an  array  of  doubles  that  will  contain  the  deviation  for 
'  setting  each  MOSZ  with  all  available  multipliers. 

t 

Sub  MinFcn(dAmbda  As  Double,  dPreset()  As  Double,  dCost()  As  Double, 
dDeviation()  As  Double) 

Dim  dCelMin  As  Double 
Dim  dCelTot  As  Double 
Dim  dCelObj  As  Double 
Dim  dCObj  As  Double 
Dim  dZJ  As  Double 
Dim  ilndex  As  Integer 
Dim  i  As  Integer 
Dim  j  As  Integer 
Dim  bContinue  As  Boolean 
dSumCost  =  0# 
dCelTot  =  0# 
dZU  =  0# 

For  i  =  1  To  (iNumMOS  -  1) 
bContinue  =  True 
Do 

If  (dPreset(i,  1)  >  -1)  Then 

dCelMin  =  dDeviation(i,  dPreset(i,  2))  +  dAmbda  *  dCost(i,  _ 
dPreset(i,  2)) 

dCObj  =  dDeviation(i,  dPreset(i,  2)) 
ilndex  =  dPreset(i,  2) 

Exit  Do 
End  If 

dCelMin  =  dlnfmity 
ilndex  =  0 

For  j  =  1  To  (iNumMult) 

If  (dDeviation(i,  j)  <=  dTest)  Then 

dZJ  =  dDeviation(i,  j)  +  dAmbda  *  dCost(i,  j) 
dCelObj  =  dDeviation(i,  j) 

If  (dZJ  <=  dCelMin)  Then 
dPreset(i,  2)  =  j 
dCelMin  =  dZJ 
dCObj  =  dCelObj 
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ilndex  =  j 
End  If 
End  If 
Next  j 

bContinue  =  False 
Loop  Until  bContinue  =  False 
dCelTot  =  dCelTot  +  dCelMin 
dZU  =  dZU  +  dCObj 

dSumCost  =  dSumCost  +  dCost(i,  ilndex) 

Next  i 

dZL  =  dCelTot  -  dAmbda  *  dBudget 
End  Sub 

'MinFeasible  is  a  subroutine  that  will  determine  minimum  feasible  object  value 
'solution  set  with  the  final  Langrangian  multiplier.  This  is  the  feasible  solution 
'set  with  the  lowest  objective  value. 

t 

'@param  -  dFLambda  is  a  double  variable  and  is  the  final  chosen  Langrangian 
multiplier. 

'  -  dSoln  is  an  array  of  doubles  that  contains  the  multipliers  solution  set. 

-  dPreset  is  an  array  of  doubles  the  list  the  preset  multipliers  if  chosen. 

-  dCost  is  an  array  of  doubles  that  will  contain  the  cost  for  setting  each 
'  MOSZ  with  all  available  multipliers. 

-  dDeviation  is  an  array  of  doubles  that  will  contain  the  deviation  for 
setting  each  MOSZ  with  all  available  multipliers. 

t 

Sub  MinFeasible(dFLambda  As  Double,  dSoln()  As  Double,  dPreset()  As  Double, 
dDeviation()  As  Double,  dCost()  As  Double) 

Dim  dCelTot  As  Double 
Dim  dCObj  As  Double 
Dim  dCelMin  As  Double 
Dim  dZJ  As  Double 
Dim  dCelObj  As  Double 
Dim  ilndex  As  Integer 
Dim  i  As  Integer 
Dim  j  As  Integer 
Dim  bContinue  As  Boolean 
dSumCost  =  0# 
dCelTot  =  0# 
dZU  =  0# 

For  i  =  1  To  (iNumMOS  -  1) 

Do 

If  (dPreset(i,  1)  >  -1)  Then 

dCObj  =  dDeviation(i,  dPreset(i,  2)) 
ilndex  =  dPreset(i,  2) 
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Exit  Do 
End  If 

dCelMin  =  din  Unity 
ilndex  =  0 

For  j  =  1  To  (iNumMult) 

If  (dDeviation(i,  j)  <=  dTest)  Then 

dZJ  =  dDeviation(i,  j)  +  dFLambda  *  dCost(i,  j) 
dCelObj  =  dDeviation(i,  j) 

If  (dZJ  <=  dCelMin)  Then 
dSoln(i)  =  j 
dCelMin  =  dZJ 
dCObj  =  dCelObj 
ilndex  =  j 
End  If 
End  If 
Next  j 

bContinue  =  False 
Loop  Until  bContinue  =  False 
dZU  =  dZU  +  dCObj 

dSumCost  =  dSumCost  +  dCost(i,  ilndex) 

Next  i 
End  Sub 

'Heuris  is  a  subroutine  that  will  improve  the  solution  to  the  Langranian  relaxation 
'of  the  SRB  multiplier  problem  by  using  a  heuristic  to  consume  any  residual  budget. 

t 

'@param  -  dDeviation  is  an  array  of  doubles  that  will  contain  the  deviation  for 
'  setting  each  MOSZ  with  all  available  multipliers. 

-  dCost  is  an  array  of  doubles  that  will  contain  the  cost  for  setting  each 
MOSZ  with  all  available  multipliers. 

-  dSoln  is  an  array  of  doubles  that  contains  the  multipliers  solution  set. 

-  dPreset  is  an  array  of  doubles  the  list  the  preset  multipliers  if  chosen. 

t 

Sub  Heuris(dDeviation()  As  Double,  dCost()  As  Double,  dSoln()  As  Double,  _ 
dPreset()  As  Double) 

Dim  dCostl  As  Double 
Dim  dCost2  As  Double 
Dim  i  As  Integer 

Dim  bRoomlnBudget  As  Boolean 
dCostl  =  dSumCost 
bRoomlnBudget  =  True 
While  bRoomlnBudget  =  True 
If  dCostl  <=  dBudget  Then 

dCost2  =NewCostl(dDeviation(),  dCost(),  dSoln(),  dPreset(),  dCostl) 

End  If 
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If  dCost2  =  dCostl  Then 
bRoomlnBudget  =  False 
Else 

dCostl  =  dCost2 
End  If 
Wend 
dZH  =  0# 
dFinalCost  =  0# 

For  i  =  1  To  (iNumMOS  -  1) 

dZH  =  dZH  +  dDeviation(i,  dSoln(i)) 
dFinalCost  =  dFinalCost  +  dCost(i,  dSoln(i)) 
dTotalDeviation  =  dTotalDeviation  +  dDeviation(i,  dSoln(i)) 

Next  i 
End  Sub 

'NewCostl  is  a  function  used  to  calculate  the  rate  of  return  for  each  MOSZ  based  on 
'the  rate  of  improvement  to  the  objective  function  that  can  be  achieve  for  each 
’dollar.  Newcostl  will  the  pick  the  MOSZ  with  the  largest  ROR  and  increase  the 
'MOSZs  solution  byone  step  in  the  multiplier  and  calculating  the  new  cost  of  the 
'solution  set. 

t 

'@param  -  dDeviation  is  an  array  of  doubles  that  will  contain  the  deviation  for 
setting  each  MOSZ  with  all  available  multipliers. 

-  dCost  is  an  array  of  doubles  that  will  contain  the  cost  for  setting  each 
MOSZ  with  all  available  multipliers. 

-  dSoln  is  an  array  of  doubles  that  contains  the  multipliers  solution  set. 

-  dPreset  is  an  array  of  doubles  the  list  the  preset  multipliers  if  chosen. 

-  dOldCostl  is  a  double  variable  that  is  the  current  cost  of  the  solution 
set. 

t 

'@return  -  NewCostl  is  a  double  variable  that  is  the  new  cost  of  the  new  solution 
set  that  has  been  heuristically  imporved. 

t 

Function  NewCostl  (dDeviation()  As  Double,  dCost()  As  Double,  dSoln()  As  Double, 
dPreset()  As  Double,  dOldCostl  As  Double)  As  Double 
Dim  dUArray()  As  Double 
Dim  dResid  As  Double 
Dim  dRNum  As  Double 
Dim  dRDen  As  Double 
Dim  dUBest  As  Double 
Dim  i  As  Integer 
Dim  j  As  Integer 
Dim  bContinuelnner  As  Boolean 
Dim  ilndex  As  Integer 
ReDim  dUArray(l  To  iNumMOS  -  1) 
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dResid  =  dBudget  -  dOldCostl 
For  i  =  1  To  (iNumMOS  -  1) 

Do  'continue  do  loop 
bContinuelnner  =  True 
If  (dPreset(i,  1)  >  -1)  Then 
dUArray(i)  =  -dlnfinity 
Exit  Do  'exits  continue  do  loop 
End  If 

If  (dSoln(i)  >=  iNuinMult)  Then 
dUArray(i)  =  -dlnfinity 
Exit  Do 

Elself  (dDeviation(i,  (dSoln(i)  +  1))  >  dTest)  Then 
dUArray(i)  =  -dlnfinity 
Exit  Do 
End  If 

dRNum  =  dDeviation(i,  dSoln(i))  -  dDeviation(i,  (dSoln(i)  +  1)) 
dRDen  =  dCost(i,  (dSoln(i)  +  1))  -  dCost(i,  dSoln(i)) 

If  ((dRDen  >  dResid)  Or  (dRNum  <=  (0.1  *  10  A  -8)))  Then 
dUArray(i)  =  -dlnfinity 
Exit  Do  ’exits  continue  do  loop 
Else 

dUArray(i)  =  dRNum  /  dRDen 
End  If 

If  (dUArray(i)  <  0#)  Then 
dUArray(i)  =  - 1  *  dUArray(i) 

End  If 

Loop  Until  bContinuelnner  =  True  'end  of  continue  loop 
Next  i 

dUBest  =  -dlnfinity 

For  i  =  1  To  (iNumMOS  -  1) 

If  (dUArray(i)  >  dUBest)  Then 
dUBest  =  dUArray(i) 
ilndex  =  i 
End  If 
Next  i 

bContinuelnner  =  True 
Do  'bContinuelnner  do  loop 
If  (dUBest  <  (-dTest))  Then 
NewCostl  =  dOldCostl 
Exit  Do  ’exits  bContinuelnner  do  loop 
End  If 

NewCostl  =  dOldCostl  -  dCost(iIndex,  dSoln(ilndex))  +  dCost(iIndex, 
(dSoln(ilndex)  +1)) 
dSoln(ilndex)  =  dSoln(ilndex)  +  1 
Loop  Until  bContinuelnner  =  True 
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End  Function 


'WriteResults  is  a  subroutine  that  will  write  the  final  solution  to  the  Opt  Solution 

'worksheet. 

t 

'@param  -  rMOSData  is  a  range  variable  that  contains  data  for  all  the  MOSZs. 

-  rBudgetlnfo  is  a  range  variable  that  contains  the  budget  and  pay 
information  needed  for  the  SRB  Optimization  problem. 

-  dSoln  is  an  array  of  doubles  that  contains  the  multipliers  solution  set. 

f 

Sub  WriteResults(rMOSData  As  Range,  rBudgetlnfo  As  Range,  dSoln()  As  Double) 
Dim  i  As  Integer 
Dim  rOutput  As  Range 

Set  rOutput  =  wsOptSoln.Range("OutputStart",  wsOptSoln.Range("OutputStart") 
.Offset(iNumMOS,  10)) 
rOutput.  Clear 
rOutput(l,  1)  =  "Occ  Field" 
rOutput(  1 ,  2)  =  "MOS" 
rOutput(l,  3)  =  "Zone" 
rOutput(l,  4)  =  "EASPop" 
rOutput(l,  5)  =  "B/Space" 
rOutput(  1,6)  =  "Multiplier" 
rOutput(l,  7)  =  "ExpReupNumber" 
rOutput(  1,8)  =  "RndReUpNumber" 
rOutput(l,  9)  =  "Over/Under" 
rOutput(l,  10)  =  "Expected  Cost" 

For  i  =  2  To  iNumMOS 

rOutput(i,  1)  =  rMOSData(i,  1) 
rOutput(i,  2)  =  rMOSData(i,  2) 
rOutput(i,  3)  =  rMOSData(i,  3) 
rOutput(i,  4)  =  rMOSData(i,  5) 
rOutput(i,  5)  =  rMOSData(i,  8) 
rOutput(i,  6)  =  (dSoln(i  -  1)  -  1)  *  dGap 

rOutput(i,  7)  =  rMOSData(i,  (rOutput(i,  6)  /  dGap  +  1 )  +  9)  *  rOutput(i,  4) 

Call  FormatTwoDecimal(rOutput(i,  7)) 

rOutput(i,  8)  =  Application.WorksheetFunction.Round(rOutput(i,  7),  0) 
rOutput(i,  9)  =  rOutput(i,  8)  -  rOutput(i,  5) 

rOutput(i,  10)  =  rBudgetInfo(3,  1)  *  rBudgetInfo((rOutput(i,  3))  +  4,  1)  *  _ 
rOutput(i,  6)  *  rOutput(i,  7) 

Call  FormatTwoDecimal(rOutput(i,  10)) 

Next  i 

End  Sub 

'File:  Utilities  Module  in  VBA 

'Date  Created:  1 7  June  2007 
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'Last  Updated:  25  July  2007 

t 

'Utilities  is  a  Module  that  contains  various  formatting  subroutines. 

t 

'@  Author  Kent  Robbins,  Jr. 

t 

Option  Explicit 

'FormatTwoDecimal  is  a  subroutine  that  will  format  a  given  range  to  two  decimal 
'places. 

f 

'@param  -  rTable  is  the  range  of  values  to  be  formatted. 

f 

Sub  FormatTwoDecimal(rTable  As  Range) 
rTable.NumberFormat  =  "0.00" 

End  Sub 

'FormatOccFieldText  is  a  subroutine  that  will  format  a  given  range  to  text  and 
'align  it  to  the  right  of  the  cell. 

f 

'@param  -  rTable  is  the  range  of  cells  to  be  formatted. 

f 

Sub  FormatOccFieldText(rTable  As  Range) 
rTable.NumberFormat  =  "@" 
rTable.HorizontalAlignment  =  xlRight 
End  Sub 

File:  WriteFiles  Module  in  VBA 
'Date  Created:  1 7  June  2007 
'Last  Updated:  25  July  2007 

t 

'WriteFiles  is  a  VBA  Modula  that  will  write  will  read  in  the  GAMS  solutions  and 
'write  the  solutions  along  with  information  about  the  MOSs  to  the  OptSoln  worksheet. 

t 

'@  Author  Kent  Robbins,  Jr. 

f 

Option  Explicit 

'WriteDataMaster  is  a  subroutine  that  will  gather  information  about  MOSs  that 
'is  needed  by  the  GAMS  solver  and  write  that  informatino  into  .dat  files. 

'@param  -  none 

t 

Sub  WriteDataMaster() 

Dim  i  As  Integer 
Dim  j  As  Integer 
Dim  rMOSData  As  Range 
Dim  rAvgPays  As  Range 
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Dim  rScalarData  As  Range 

Dim  sFileNames()  As  String 

Dim  slndex()  As  String 

Dim  iLength  As  Integer 

Dim  dMaxTrnCost  As  Double 

ChDrive  Left$(ActiveWorkbook.Path,  1) 

ChDir  ActiveWorkbook.Path 

Set  rMOSData  =  wsReducedProblem.Range(wsReducedProblem.Range  _ 
("RFPOutputStart").End(xlToRight),  wsReducedProblem.Range  _ 
("RFPOutputStart").End(xlDown)) 

Set  rAvgPays  =  wsReducedProblem.Range("AvgPayA",  wsReducedProblem  _ 
.Range("AvgPayA").Offset(3,  0)) 

Set  rScalarData  =  wsReducedProblem.Range("OverageWeight",  wsReducedProblem 
.Range("OverageWeight").Offset(3,  0)) 
iLength  =  rMOSData.Rows. Count 
dMaxTrnCost  =  0# 

ReDim  sFileNames(l  To  15)  As  String 
ReDim  slndex(l  To  iLength)  As  String 
sFileNames(l)  =  "MOSZoneList.dat" 
sFileNames(2)  =  "MultiplierList.dat" 
sFileNames(3)  =  "AvgPay.dat" 
sFileNames(4)  =  "MOSPop.dat" 
sFileNames(5)  =  "EASPop.dat" 
sFileNames(6)  =  "PresetMult.dat" 
sFileNames(7)  =  "WgtFactor.dat" 
sFileNames(8)  =  "BoatSpace.dat" 
sFileNames(9)  =  "TrnCost.dat" 
sFileNames(lO)  =  "ReenlistmentRates.dat" 
sFileNames(l  1)  =  "MultiplierNumbers.dat" 
sFileNames(12)  =  "Q.dat" 
sFileNames(13)  =  "Budget.dat" 
sFileNames(14)  =  "Yrs.dat" 
sFileNames(15)  =  "MaxTrnCost.dat" 

For  i  =  2  To  iLength 

If  (rMOSData(i,  3)  =  1)  Then 

slndex(i)  = . +  CStr(rMOSData(i,  2). Value)  +  "a"  + . 

Elself  (rMOSData(i,  3)  =  2)  Then 

slndex(i)  = . +  CStr(rMOSData(i,  2). Value)  +  "b"  + . 

Else 

slndex(i)  = . +  CStr(rMOSData(i,  2). Value)  +  "c"  + . 

End  If 

If  (rMOSData(i,  9)  >  dMaxTrnCost)  Then 
dMaxTrnCost  =  rMOSData(i,  9) 

End  If 
Next  i 
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For  i  =  1  To  1 1 

Call  WriteData(i,  sFileNames(i),  slndex(),  rMOSData,  rAvgPays,  iLength, 
rScalarData) 

Next  i 

For  i=  12  To  15 

Dim  dSend  As  Double 
If  i  =  15  Then 

dSend  =  dMaxTrnCost 
Else 

dSend  =  rScalarData(i  -  1 1) 

End  If 

Call  WriteDataScalar(sFileNames(i),  dSend) 

Next  i 
End  Sub 

'WriteDataScalar  is  a  subroutine  that  will  write  scalar  data  into  a  .dat  file  for 
’GAMS. 

t 

'@param  -  none 

t 

Sub  WriteDataScalar(sNames  As  String,  dData  As  Double) 

Dim  iFN  As  Integer 
iFN  =  FreeFile() 

Open  sNames  For  Output  As  iFN 
Print  #iFN,  7"  &  dData  &  7" 

Close  iFN 
End  Sub 

'WriteData  is  a  subroutine  that  will  write  MOSZ  data  into  a  .dat  file  for 
'GAMS. 

t 

'@param  -  iPosition  As  Integer  is  the  index  for  the  files  to  be  written. 

-  sNames  As  String  is  the  file  name  for  the  data  to  be  written. 

-  slndexQ  As  String  is  the  name  of  the  MOSZ  to  be  written. 

-  rMOSData  As  Range  is  the  range  of  data  to  be  printed. 

'  -  rPay  As  Range  variable  of  the  avg  pays  to  be  prin  ted. 

-  iLength  As  Integer  is  the  number  of  files  to  be  printed. 

-  rScalarData  As  Range  variable  of  all  the  scalar  values  needed. 

t 

Sub  WriteData(iPosition  As  Integer,  sNames  As  String,  slndex()  As  String,  _ 
rMOSData  As  Range,  rPay  As  Range,  iLength  As  Integer,  _ 
rScalarData  As  Range) 

Dim  i  As  Integer,  j  As  Integer,  iFN  As  Integer,  iNumMult  As  Integer 
iFN  =  FreeFile() 

Open  sNames  For  Output  As  iFN 
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If  (iPosition  =1)  Then 
Print  #iFN,  7" 

For  i  =  2  To  iLength 
Print  #iFN,  slndex(i); 

Ifi  <>  iLength  Then 
Print  #iFN, 

End  If 
Next  i 

Print  #iFN,  7" 

Elself  (iPosition  =  2)  Then 
Print  #iFN,  7" 

For  i  =  1  To  (rMOSData.Columns. Count  -  (9)) 
Print  #iFN,  i; 

If  i  <>  (rMOSData.Columns. Count  -  (9))  Then 
Print  #iFN, 

End  If 
Next  i 

Print  #iFN,  7" 

Elself  (iPosition  =  3)  Then 
Print  #iFN,  7" 

For  i  =  2  To  iLength 
Print  #iFN,  slndex(i)  & 

If  (rMOSData(i,  3)  =  1)  Then 
Print  #iFN,  rPay(l) 

Elself  (rMOSData(i,  3)  =  2)  Then 
Print  #iFN,  rPay(2) 

Else 

Print  #iFN,  rPay(3) 

End  If 

If  (i  <>  iLength)  Then 
Print  #iFN, 

End  If 
Next  i 

Print  #iFN,  7" 

Elself  (iPosition  =10)  Then 

iNumMult  =  rMOSData.Columns. Count  -  9 
Print  #iFN,  "dummy"; 

For  i  =  1  To  iNumMult 
Print  #iFN,  &  i; 

Next  i 
Print  #iFN, 

For  i  =  2  To  iLength 
Print  #iFN,  slndex(i); 

For  j  =  iPosition  To  rMOSData.Columns. Count 
Print  #iFN,  &  rMOSData(i,  j); 
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Next  j 
Print  #iFN, 

Next  i 

Elself  (iPosition  =11)  Then 

iNumMult  =  rMOSData. Columns. Count  -  9 
Print  #iFN,  "/" 

For  i  =  1  To  iNumMult 

Print  #iFN,  i  &  "="  &  rMOSData(l,  9  +  i); 

If  (i  <>  iNumMult)  Then 
Print  #iFN, 

End  If 
Next  i 

Print  #iFN,  "/" 

Elself  (iPosition  =12)  Then 
Print  #iFN,  "/" 

For  i  =  1  To  4 

Print  #iFN,  i  &  "="  &  rScalarData(i); 

If  (i  <>  4)  Then 
Print  #iFN, 

End  If 
Next  i 

Print  #iFN,  7" 

Else 

Print  #iFN,  "/" 

For  i  =  2  To  iLength 

Print  #iFN,  slndex(i)  &  "="  &  rMOSData(i,  iPosition); 

If  (i  <>  iLength)  Then 
Print  #iFN, 

End  If 
Next  i 

Print  #iFN,  "/" 

End  If 
Close  iFN 
End  Sub 

'File:  Solve  Module  in  VBA 
'Date  Created:  1 7  June  2007 
'Last  Updated:  30  July  2007 

t 

'SolveGAMS  is  a  subroutine  that  will  run  the  GAMS  solution  model  in  a  shell 
'program. 


'@  Author  Kent  Robbins,  Jr. 

t 

Option  Explicit 
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Sub  SolveGAMS() 

ChDrive  Left$(ActiveWorkbook.Path,  1) 

ChDir  ActiveWorkbook.Path 

Shell  "gams  ThesisSRBModel.gms",  vbNormalFocus 
End  Sub 

'File:  ReadResuIts  Module  in  VBA 
'Date  Created:  1 7  June  2007 
'Last  Updated:  14  July  2007 

f 

'ReadResuIts  is  a  VBA  Module  that  will  read  in  the  optimum  SRB  multiplier  results  from 
'GAMS  and  print  them  to  Opt  Solution  worksheet  along  with  its  associated  data. 

t 

'@  Author  Kent  Robbins,  Jr. 

t 

Option  Explicit 

'GAMSResults  is  a  subroutine  that  will  read  in  the  optimum  SRB  multiplier  results  as 
'an  outputfrom  GAMS.  Results  will  be  outputted  with  other  important  data  for  the 
'MOS's. 

t 

'@param  -  None 

t 

Sub  GAMSResults() 

Dim  iFN  As  Integer 

Dim  i  As  Integer 

Dim  j  As  Integer 

Dim  iLength  As  Integer 

Dim  r  As  Double 

Dim  dGap  As  Double 

Dim  dExpReUpNum  As  Double 

Dim  rMOSData  As  Range 

Dim  rOutput  As  Range 

Dim  rAvgPays  As  Range 

ChDrive  Left$(ActiveWorkbook.Path,  1) 

ChDir  ActiveWorkbook.Path 

Set  rMOSData  =  wsReducedProblem.Range(wsReducedProblem.Range  _ 
("RFPOutputStart").End(xlToRight),  wsReducedProblem.Range  _ 
("RFPOutputStart").End(xlDown)) 

Set  rAvgPays  =  wsReducedProblem.Range("AvgPayA",  wsReducedProblem  _ 
.Range("AvgPayA").Offset(3,  0)) 
iLength  =  rMOSData.Rows. Count 
dGap  =  rMOSData(l,  1 1)  -  rMOSData(l,  10) 

Set  rOutput  =  wsOptSoln.Range("OutputStart",  wsOptSoln.Range("OutputStart")  _ 
.Offset(  iLength,  10)) 
rOutput(l,  1)  =  "Occ  Field" 
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rOutput(  1 ,  2)  =  "MOS" 
rOutput(l,  3)  =  "Zone" 
rOutput(l,  4)  =  "EASPop" 
rOutput(l,  5)  =  "B/Space" 
rOutput(  1,6)  =  "Multiplier" 
rOutput(l,  7)  =  "ExpReupNumber" 
rOutput(l,  8)  =  "RndReUpNumber" 
rOutput(l,  9)  =  "Short/Over" 
rOutput(l,  10)  =  "Exact  Expected  Cost" 
iFN  =  FreeFile() 

Open  "SRBOptResults.csv"  For  Input  As  iFN 
For  i  =  2  To  iFength 

rOutput(i,  1)  =  rMOSData(i,  1) 
rOutput(i,  2)  =  rMOSData(i,  2) 
rOutput(i,  3)  =  rMOSData(i,  3) 
rOutput(i,  4)  =  rMOSData(i,  5) 
rOutput(i,  5)  =  rMOSData(i,  8) 

Input  #iFN,  r 

dExpReUpNum  =  rMOSData(i,  (r  /  dGap  +  1)  +  9)  *  rOutput(i,  4) 

If  dExpReUpNum  =  0  Then 
rOutput(i,  6)  =  0 

dExpReUpNum  =  rMOSData(i,  10)  *  rOutput(i,  4) 

Else 

rOutput(i,  6)  =  r 
End  If 

rOutput(i,  7)  =  dExpReUpNum 
Call  FormatTwoDecimal(rOutput(i,  7)) 

rOutput(i,  8)  =  Application.WorksheetFunction.Round(rOutput(i,  7),  0) 
rOutput(i,  9)  =  rOutput(i,  8)  -  rOutput(i,  5) 

rOutput(i,  10)  =  4  *  rAvgPays(rOutput(i,  3))  *  rOutput(i,  6)  *  rOutput(i,  7) 
Call  FormatTwoDecimal(rOutput(i,  10)) 

Input  #iFN,  r 
Next  i 
Close  iFN 
End  Sub 
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APPENDIX  C:  EXAMPLE  INPUTS  AND  OUTPUTS 


This  appendix  contains  the  examples  of  the  inputs  and  outputs  of  the  Excel 
interface  that  was  used  to  solve  the  SRB  multiplier  selection  problem  for  FY04. 
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Figure  1 .  Excel  Worksheet  for  Zone  A  Reenlistment  Data. 
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Figure  2.  Excel  Worksheet  of  complied  MOSZs  and  data 
need  for  the  execution  of  the  Model. 
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Figure  3.  Excel  Worksheet  with  MOSZ  data,  multiplier 

solution  and  expected  cost. 
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